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Abstract 

Searches were performed for topologies predicted by gauge-mediated Supersymmetry break- 
ing models (GMSB). All possible lifetimes of the ncxt-to-lightest SUSY particle (NLSP), either 
the lightest ncutralino or slcpton, decaying into the lightest SUSY particle, the gravitino, were 
considered. No evidence for GMSB signatures was found in the OPAL data sample collected at 
centre-of-mass energies of -y/s — 189-209 GeV at LEP. Limits on the product of the production 
cross-sections and branching fractions are presented for all search topologies. To test the impact 
of the searches, a complete scan over the parameters of the minimal model of GMSB was per- 
formed. NLSP masses below 53.5 GeV/c^ in the neutralino NLSP scenario, below 87.4GeV/c^ 
in the stau NLSP scenario and below 91.9 GeV/c^ in the slepton co-NLSP scenario are excluded 
at 95 % confidence level for all NLSP lifetimes. The scan determines constraints on the universal 
SUSY mass scale A from the direct SUSY particle searches of A > 40, 27, 21, 17, 15TeY/c^ for 
messenger indices A = 1, 2, 3, 4, 5 for all NLSP lifetimes. 
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1 Introduction 



Supersymmetry (SUSY) [1] provides a method of solving the hierarchy problem by introducing a 
set of new particles which cancel the large radiative corrections to the Higgs mass. The cancellation 
is achieved by assuming that, for each Standard Model (SM) particle chirality state, there is one 
additional particle identical to its SM partner except that its spin differs by 1/2 unit. If SUSY were 
an exact symmetry, the new SUSY particles would have the same masses as their SM partners. 
Since this scenario is experimentally excluded, SUSY must be a broken symmetry. It is typically 
assumed that SUSY is broken in some "hidden" sector of new particles, and is "communicated" 
(or mediated) to the "visible" sector of SM and SUSY particles by one of the known interactions. 
The two scenarios for this mediation that have been most widely investigated are gravity and gauge 
mediation. 

In gauge-mediated SUSY breaking (GMSB), the hidden sector can lie at masses as low as about 
10^ GeV/c^. In most current GMSB theoretical work [2-4], it is assumed that this sector is cou- 
pled to a messenger sector, which in turn couples to the visible sector through normal SM gauge 
interactions. In its minimal version five new parameters and a sign are introduced in addition to 
the SM parameters, usually chosen to be the SUSY breaking scale the messenger scale M, the 
messenger index N, the ratio of the vacuum expectation values of the two Higgs doublets tan/3, 
the sign of the Higgs sector mixing parameter sign(//), and the mass scale A, which determines the 
SUSY particle masses at the messenger scale. 

A feature which distinguishes gravity-mediated from gauge-mediated SUSY breaking models is 
the mass of the gravitino, G. In gravity-mediated models, the G is usually too heavy to have a signif- 
icant effect on SUSY phenomenology, while in GMSB models the G is typically light (< 1 MeV/c^) 
and is the lightest SUSY particle, the LSP. While the G is a spin 3/2 particle, only its it 1/2 spin 
projections (which have absorbed the goldstino associated with spontaneous SUSY breaking via the 
"superhiggs" mechanism [5]) interact with weak, rather than gravitational, strength interactions, 
and contribute to phenomenology. 

In GMSB, the next-to-lightest SUSY particle (NLSP) is either the lightest neutralino, Xi (neu- 
tralino NLSP scenario), or the lightest scalar lepton^, 1^, which is either the mass degenerate e/j, 
fLf{ and fi (slepton co-NLSP scenario) or in a significant fraction of the parameter space the lightest 
scalar tau lepton, fi (stau NLSP scenario). For SUSY particles heavier than the NLSP, the coupling 
to the G is very small, and typically they will decay to the NLSP, which then decays to the gravitino 
via xi ^ tG or ^ ^ iG. 

One feature of GMSB is that the decay length f3^cT of the NLSP depends on the intrinsic SUSY 
breaking scale ^/F : 

\ ifiNLSP J l^lOOTeVy \m%^sp J 

where mNLSP and -Enlsp are mass and energy of the NLSP. Due to the range allowed for ^/F, from 
2000TeV/c^ down to lOOTeV/c^ [6], the decay length is effectively unconstrained and decays inside 
and outside the detector volume have to be considered. If the NLSP decay to the gravitino occurs 
with a small lifetime, event signatures will include energetic leptons or photons, plus significant 

^For the first and second generations in tlie minimal GMSB model, the lighter state is predominantly right-handed, 
since the off-diagonal components of the mass matrix are small due to small Yukawa couplings. Therefore in the 
following for the selectron and the smuon the symbol £r instead of £i is often used. 
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missing energy, ^, due to the undetected gravitinos. For intermediate lifetimes, the events may 
include tracks not pointing to the primary interaction point (tracks with large impact parameters) 
or kinked tracks. For long lifetimes, the signatures can include tracks from heavy long-lived charged 
particles. 

In this paper, we present results on searches for GMSB signatures with either a or a £i 
(or fi) NLSP using the data sample acquired by the OPAL detector at centre-of-mass energies of 
^/s = 189-209 GeV. All relevant signatures - direct NLSP pair-production and its appearance in the 
decay chains of heavier SUSY particles like charginos, neutralinos or sleptons - for all NLSP lifetimes 
are considered in our searches. Topologies and sparticle decay modes together with a reference to 
the description of the analyses are given in Table 1 for the slepton and stau NLSP scenarios and in 
Table 2 for the neutralino NLSP scenario. 

The paper is organized as follows. The OPAL detector, data and Monte Carlo samples are 
introduced in Section 2, followed by a brief description of the event selections of different topologies 
in the slepton and stau NLSP scenarios in Section 3 and for the neutralino NLSP scenario in 
Section 4. The strategy to combine such a large number of analyses is described in Section 5. 
Results are presented in Section 6, and finally in Section 7 constraints on the parameters of the 
minimal GMSB model are discussed. 

2 Data and Monte Carlo samples 

A detailed description of the OPAL detector can be found elsewhere [12]. Tracking of charged 
particles was performed by a central detector, enclosed in a solenoid which provided a uniform 
axial magnetic field of 0.435 T. The central detector consisted of a two-layer silicon microvertex 
detector [13], a high precision vertex chamber with both axial and stereo wire layers, a large volume 
jet chamber providing both tracking and ionization energy loss information, and additional chambers 
to measure the z coordinate of tracks as they leave the central detector^ . These detectors provided 
tracking coverage for polar angles | cos 6\ < 0.96, with a typical transverse momentum (pt) resolution 
of (7p,j,/p^ = 1.25 X 10^'^ {GeV / c)~^ . The solenoid coil was surrounded by a time-of-flight counter 
array and a barrel lead-glass electromagnetic calorimeter with a pre-sampler. Together with the 
endcap electromagnetic calorimeters, the lead-glass blocks covered the range | cos 6*1 < 0.98. Outside 
the electromagnetic calorimetry, the magnet return yoke was instrumented with streamer tubes 
to form a hadronic calorimeter with angular coverage in the range |cos^| < 0.91. The region 
0.91 < [cos^l < 0.99 was instrumented with an additional pole-tip hadronic calorimeter using 
multi-wire chambers. The detector was completed with muon detectors outside the magnet return 
yoke. These were composed of drift chambers in the barrel region and limited streamer tubes in the 
endcaps, and together covered 93 % of the full solid angle. 

The data used for the new particle searches described here were taken at centre-of-mass energies 
between 189 GeV and 209 GeV during the LEP running periods from 1998 to 2000. The total inte- 
grated luminosity is about 600 pb~^, evaluated using small angle Bhabha scattering events observed 
in the forward calorimeters [14] with an error of about 0.3%. The data were analyzed in 8 inde- 
pendent bins in i/s, summarized in Table 3. Large samples of Monte Carlo simulated events were 

right-handed coordinate system is adopted, where the a;-axis points toward the centre of the LEP ring and the 
z-axis points in the direction of the electron beam. The polar angle 6 and the azimuthal angle (p are defined w.r.t. z 
and X, respectively. 
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Table 1: GMSB signatures for the stau NLSP and the slepton co-NLSP scenario. The notation (1) 
"prompt", (2) "within detector" and (3) "outside detector" refers to a NLSP decay such that the 
decay vertex is (1) close to the interaction region and not measurably displaced, (2) resolvable by 
the large impact parameter and kinked track searches and (3) well outside the detector. For each 
signature a reference to the description of the analysis is given, either to a section of this or to a 
separate paper. 
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Table 2: GMSB signatures for the neutralino NLSP scenario. The notation (1) "within detector" 
and (2) "outside detector" refers to a NLSP decay such that the decay vertex is (1) at the interaction 
region (prompt decay) up to a few tenths of centimetres inside the detector volume and (2) well 
outside. For each signature a reference to the description of the analysis is given, either to a section 
of this or to a separate paper. 

generated at each ^/s to optimize the search algorithms for new particles and were used to evaluate 
their efficiency and the number of expected events from Standard Model sources. 

Signal events were generated for all search topologies with 10 different NLSP lifetimes ranging 
from a very short lifetime, r = 10~^^ s, corresponding to a prompt NLSP decay with a decay length 
of 300 ^m for a particle with (5^ = p/m = 1, up to very long lifetimes, r = 10~^ s, corresponding to 
a stable NLSP with a decay length of 300 m which always decays outside the detector. Mass grids 
from 45 GeV/c^ to the kinematic limit with 10-100 points per channel per y/s bin and 1000 events 
per grid point were used for the full simulation of signal Monte Carlo events. The SUSYGEN [15] 
event generator was used to simulate the signal events for and Xi pair-production. For chargino 
pair-production the W and Z boson widths can play an important role and are not fully treated 
in SUSYGEN. The DFGT generator [16] was used to simulate these signal events. It includes spin 
correlations and allows for a proper treatment of the W boson and the Z boson width effects in 
the chargino decay. Both SUSYGEN and DFGT include initial-state radiation (ISR). The Lund 
string fragmentation scheme of PYTHIA in the JETSET package [17] was used to describe the 
hadronization. The gravitino mass was set to zero in the generation, since a small mass up to 
©(MeV) favoured by the models has a negligible effect on the detection efficiencies. The GEANT3 
[18] software package, which was used to simulate the transport of particles through the material 
of the detector, only allows isotropic two- or three-body decays of long-lived particles. Thus an 
interface between GEANT and the event generators was created to obtain correct angular and 
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energy distributions for the decay products of the NLSP, which decays in flight. 

The simulated Standard Model background is composed of two-photon (e"'^e~ e~''e~77 — > 
e^e~ ff), four-fermion (e'''e~ — > ////), multi-photon (e^e~ 77(7)) and two-fermion (e~'"e~ — > 
(Z/7)* II and e+e" — > (Z/7)* qq) final states. Leptonic two-photon processes were simulated 
with the BDK [19] and Vermaseren [20] programs, and hadronic two-photon processes with the 
PHOJET [21] and HERWIG [22] programs. To simulate four-fermion processes, the KORALW [23] 
generator (which simulates ISR using a spectrum with transverse momentum) and, for final states 
with electrons, the grc4f [24] program (with collinear ISR) were used. Multi-photon final states were 
simulated with the RADCOR [25] program. For Bhabha events, the BHWIDE [26] and TEEGG [27] 
generators were used. Muon- and tau-pairs were simulated using the KK2F [28] program and I'l'i^) 
events with the NUNUGPV [29] generator. Multi-hadronic events were simulated using PYTHIA 
and KK2F. The equivalent luminosities of the background Monte Carlo samples were in most cases 
much larger than the luminosity of the data; however, in a few background samples (especially 
the high cross-section two-photon processes), the equivalent luminosity of the background Monte 
Carlo was only about 2x that of the data. Final states with six or more primary fermions were 
not included in the background Monte Carlo samples but they are expected to make a negligible 
contribution. 

The simulated signal and background events were processed through the full simulation of the 
detector, and the same analysis chains were applied to Monte Carlo events as to the data. The 
statistical errors on the signal and background Monte Carlo are included in all results. All search 
analyses described in this paper used the same sets of simulated background events. 



3 Searches in the stau NLSP and slepton co-NLSP scenarios 

In the scenarios where mass-degenerate charged sleptons or the stau are the next-to-lightest SUSY 
particles, they decay to their charged Standard Model partner and the electrically neutral, lightest 
SUSY particle, the gravitino. Depending on the NLSP lifetime, different analyses were applied. In 
the following section all searches for SUSY particles in the stau NLSP and the slepton co-NLSP 
scenario are presented. In Table 3 all search results not described in other papers are summarized. 

3.1 Searches for promptly decaying charged NLSPs 

The decay vertex of an NLSP with a lifetime less than 10^^^ s is close to the interaction region and 
considered not to be displaced from the primary vertex. Depending on whether the NLSP is directly 
pair-produced or appears in the decay chain of heavier SUSY particles, two or more charged leptons 
are expected in an event and different search algorithms were optimized to select candidate events 
with different track multiplicities. 

3.1.1 Search for di-lepton events with missing transverse momentum 

Candidate events for pair-produced NLSPs, decaying promptly into a lepton and a gravitino (e+e~ — > 
^R^lii '^R ^^), should have the signature of two leptons plus missing transverse momentum. The 
search for events of this type is described in [7]. The number of di-lepton events, the dependence 
on centre-of-mass energy and the event properties are consistent with expectations from Standard 
Model processes, predominantly W+W" production with both W bosons decaying leptonically. 
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No evidence for new phenomena is apparent. The topology of two leptons plus missing transverse 
momentum can also result from the pair-production of charginos which decay to an invisible neutrino 
and the charged NLSP (e+e" — > XiXi ^ xf ~^ ^R^^ ^G). Thus the same analysis as for pair- 

produced NLSPs is applied. 

3.1.2 Selection of events with more than two leptons and missing energy 

The NLSP might appear in the decay chain of a heavier pair-produced SUSY particle, a neutralino 
or a heavy slepton. In such cases in addition to the two leptons from the prompt decays of both 
NLSPs, two or four more leptons are expected in the detector (e'^e" XiXii Xi ~^ ^R^^ ^R ~^ 
or e+e- ^ £+1]^, Ir ^ ^x?, X? ^ rf, f ^ tG). 

An analysis to search for such topologies with multi-leptons and missing energy which performs 
no explicit electron or muon identification was designed to maintain simplicity, minimize systematic 
errors and maximize efficiency for all lepton flavours. After the event reconstruction, double-counting 
of energy between tracks and calorimeter clusters was corrected by reducing the calorimeter cluster 
energy by the expected energy deposition from aligned charged tracks [8], including particle identi- 
fication information. Charged tracks and neutral clusters were defined to be of "good" quality using 
the requirements of [8]. 

The analysis was split into three steps: (1) preselection of events with low to intermediate 
multiplicity and missing energy; (2) additional kinematic requirements to remove unmodeled back- 
grounds; (3) requirements that each event contains at least four jets, each consistent with being an 
isolated lepton (including a tau). 

(1) Preselection 

(PI) The number of good charged tracks A'^trk in the event satisfied 4 < A'^ti-k < 10. This cut removed 
a large fraction of the lepton pair and multihadronic backgrounds. 

(P2) Events that may include charged tracks originating from noise in the central jet chamber were 
removed by vetoing events containing any tracks with dE'/dx < 5 keV/cm. 

(P3) A'trk/A^triJ > 0.20, where N'^l^ is the total number of charged tracks reconstructed in the event. 
This reduced the number of events due to machine backgrounds (eg', beam-gas and beam- wall 
interactions). 

(P4) Cosmic rays were vetoed using information from the tracking and time-of-flight systems. 

(P5) Mvis > 5.0 GeV/c^ where Mvis is the visible mass of the event. This cut was required to 
remove events which were not simulated in the untagged two-photon Monte Carlo generation. 

(P6) Pt. miss > 0.05 X -y/s, where -PT,miss is the magnitude of the component of event momentum 
transverse to the beam axis. This greatly reduced the number of two-photon and two-fermion 
events. 

(P7) Cuts were made on the maximum energy deposited in either side of the forward calorimeters 
in order to reduce two-photon and some machine backgrounds. These cuts have an associ- 
ated efficiency loss due to activity which may be present in these subdetectors for signal (or 
background) events, which was evaluated with random beam-crossing triggers. The typical 
efficiency loss was about 3%. 
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(2) Kinematic requirements 

There is an excess of events surviving the preselection. This excess does not have the characteristics 
of any of the signal hypotheses considered in this paper, and it was removed with the following 
requirements. 

(Kl) (pacop > 10°. The acoplanarity angle (/^acop was calculated by forcing the event into two jets 
using the Durham jet finder [30], and subtracting the opening angle between the two jets in 
the plane transverse to the beam axis from 180°. This cut removed some surviving background 
from lepton pair production and two photon events. 

(K2) I cos^missl < 0.90, where ^miss is the polar angle of the event missing momentum vector. This 
cut removed some residual surviving Bhabha scattering and two photon background. 

(K3) 0.20 < Eyis/^/s < 0.90 where E'vis is the event visible energy. This cut primarily removed 
two-photon background. The distributions of the event visible energy in data and simulated 
background are illustrated in Figure 1 (a) and (b) . 

(3) Lepton-jet identification 

The events surviving the preselection and kinematic requirements were split up into jets. The jets 
were then required to be consistent with having originated from electrons, muons, or tau-leptons. 
The signal hypotheses considered here are most often dominated by tau-leptons, and no explicit 
electron or muon identification tools are used. The event was split into jets using the Durham jet 
finder [30] with ycut = 0.0004. This ycut was chosen as the smallest value (to maximize efficiency for 
finding leptons which are close to each other) which would not split single 3-prong tau-lepton jets 
into multiple jets. The following cuts were then applied. 

(1) iVf-jct ^ 4, where Ni_j(,t is the number of jets in the event compatible with being from a 
parton-level lepton: 

- Lepton-jet momentum satisfies 0.01 < P£_jct/^bcam < 0.80; 

- Lepton-jet mass satisfies Mg.jct < 5 GeV/c^; 

- Number of charged tracks satisfies 1 < A'trk ^ 3. 

The distributions of the number of lepton-jets in data and simulated background are illustrated 
in Figure 1 (c) and (d). 

(2) ^isoi > 10° where ^isoi is the smallest opening angle between any two lepton-jets. 

(3) S| cos0£_jet| < 3, where S| cos0£_jct| is the sum of the magnitudes of the cosines of the polar 
angles of the lepton-jet candidates. This cut removed a small residual amount of e^e~ — > qq 
background. 

(4) E'sumlnon £ — jet)/Evis < 0.20, where E'sumlnon ^ — jet)/Evis is the fraction of the event visible 
energy which is not contained in a lepton-jet candidate. 

The full analysis retains efficiency for the GMSB processes XiXi ■, both for the case where 

the neutralino decays with equal branching ratios into all three slepton generations (equal BR), and 
also when it decays to fr with a 100 % branching fraction. Typically when the mass difference 
M(xi) — M(€) is above 5GeV/c^, the efficiency for the equal BR scenario is 40-50%, while for the 
100% fr decay it is 30-40%. For M(xi) — M{1) = 5 GeV/c^, the corresponding efficiencies are 
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typically 30% and 15%. For M{xi) — M{£) as low as 2 GeV/c^, the equal BR scenario retains 
about 10% efficiency, while the very soft recoil tau leptons make the 100% fr BR have about 1% 
efficiency. The full selection yields two events in the data with 2.8 it 0.3(stat)^Q3(syst) expected 
from the background. 

To search for selectron and smuon pair-production, the final analysis cut (on the fraction of 
energy in the event which is not included in a lepton candidate jet) has relatively poor efficiency. 
This is due to the 6 leptons in the event, any of which could fail the lepton-jet classification. The 
optimal search for that process therefore excludes the final cut, and selects 5 events in the data 
with 4.2 =b 0.3(stat)^Q 4(syst) expected in the background. The mass grid in this mode is three 
dimensional, since the decay proceeds via 1 £x^ — > ifr, and therefore the neutralino mass is also 
important with the efficiency depending strongly on both M{i) — M(xi) and M{xi) ~ M{f). The 
detection efficiency for e'''e~ and are similar, with the /U"*"/!" process always slightly higher. 

If all mass differences in the process are >5 GeV/cP, the detection efficiencies are about 50%. For 
M{i) — M{xi) as low as 2 GeV/c^ (soft electrons or muons) the efficiency is about 35-45%, while 
for M(xi) — M(f) as low as 2 GeV/c^ (soft taus) the efficiency is typically reduced to about 15%. 

3.1.2.1 Systematic uncertainties 

The dominant background in the selection is from non-multiperipheral 4-fermion production 
(which in our simulations is almost entirely due to the production of four charged leptons), consti- 
tuting 80 % of the total background expectation. A smaller contribution of about 15 % is expected 
from the multiperipheral diagrams e^e~ — > e^e"^"*"^" (about 95% from e~''e~ — > e~^e~T'^T~). The 
remaining 5% is from e^e~ t^t~ . All other background processes are negligible. The Monte 
Carlo modeling of the background was tested using different event generators. Of particular concern 
is the fact that, in our principal background samples available at all centre-of-mass energies, the 
multiperipheral diagrams are treated with specialized "two-photon" event generators, which neglect 
interference with non-multiperipheral diagrams. Special samples of the full set of e"'"e~ — > e~^e~i~^£~ 
diagrams, including interference, were prepared using grc4f2.2 [24] at ^/s = 206 GeV to study this 
effect. The background using the full set of e'^e~i~^i~ diagrams including interference changes by 
(—41 lb 5) % after the preselection compared to using our standard set of Monte Carlo generators. 
After all cuts, the result is (-1-20 it 20) %. The systematic error for this effect is taken to be 40% on 
the e'^e~i~^£~ component of the background. The modeling of the 4-fermion background, excluding 
e~^e~i~^i~ , was tested by comparing grc4f to our primary event generator K0RALW4f. After the 
preselection, the background from grc4f changes by (-1-11 it 4) %, while after all cuts it changes by 
(—10 lb 30) %. A 10% systematic error is assigned to the non e~^e~i~^l~ 4-fermion background from 
this check. The modeling of e"'"e~ t~^t~ was checked comparing the predictions of KK2f and 
KoralZ. The difference was negligible (< 1 %) and is neglected. 

The modeling of the cut variables by the detector simulation can introduce additional system- 
atic errors. This modeling was checked by comparing data samples with "relaxed cuts", in which 
any signal contribution is expected to be negligible, with Monte Carlo simulation to determine the 
uncertainty on each cut variable. Each cut position was varied in each direction by this uncer- 
tainty, and the background level and signal efficiency re-evaluated. The lepton-jet identification was 
checked using different jet finders and jet resolution parameters, but there was no indication for 
any systematics from that procedure and no additional errors are assigned. From these studies, the 
systematic error on the background due to the Monte Carlo modeling of the cut variables is taken as 
-I-31.5/-5.6 %, and on the relative signal efficiency as -I-4.2/-4.0 %, which is symmetrised to ±4%. 
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3.2 Searches for charged NLSPs with short Hfetime (tracks with large impact 
parameters) 



If the NLSP is a slepton and its lifetime is such that the decay happens within a few centimetres 
of the primary interaction point, no NLSP track can be reconstructed. Only secondary tracks of 
its decay products are visible and they have large impact parameters with respect to the primary 
interaction point. In the following section the analyses searching for large impact parameter tracks 
consistent with short-lived NLSPs with lifetimes between roug hly 10-1^ s and lO'^ s are discussed. 

All search topologies expected from SUSY particles in the stau or slepton co-NLSP scenarios 
have in common that, for short NLSP lifetimes, at least two tracks with large impact parameters 
would be present. The number of additional tracks depends on the pair-produced SUSY particle. 
The number of good tracks [8] in the event was determined, and then all tracks in the event were 
classified into one of three categories. 

Primary tracks: Good tracks with a distance from the point of closest approach to the event 
vertex in the x-y plane, the impact parameter |(io|, smaller than 0.05 cm. 

Secondary tracks: Good tracks with an impact parameter |(io| greater than 0.05cm, a high 
impact parameter precision of (5do/|c^o| < 0.1 and a significant amount of transverse momentum 
(PT > 1.5GeV/c). 

Additional tracks: All remaining tracks. 
Candidate events must contain at most ten good tracks and at least two secondary tracks. 

3.2.1 Event selection for direct NLSP and chargino pair-production 

For a smuon or selectron NLSP, the processes e+e~ — > l^i^ and e+e" xtxi (xf ~^ ^K^) followed 
by the delayed decay of the NLSP (^r —>■ iG) would lead to a signature of two tracks in the event, 
both of which can have large impact parameters. 

(1) It was required that exactly two secondary tracks, no primary tracks and no more than three 
additional tracks were found. 

After this topological cut, a selection was applied to separate the signal from the background, 
mainly two-photon or (radiative) Bhabha events, cosmic particles traversing the detector, hadronic 
interactions with the detector material or beam wall and beam gas events. 

(2) To reduce combinatorial background, the two secondary tracks were required to have different 
charges to satisfy charge conservation. 

(3) The barrel time-of-flight detector and the central jet chamber were used to reject cosmic ray 
particles. 

(4) To reduce the background from beam-gas events, which often contain tracks coming from an 
interaction point that is shifted by a large amount along the z-axis with respect to the main 
event vertex, the longitudinal impact parameter zq of each secondary track was required to 
satisfy \zq\ < 40cm. If both tracks had a \zq\ greater than 6cm, the longitudinal impact 
parameters zq of the two tracks were required to have opposite signs. 
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(5) To reject background containing back-to-back tracks, especially Bhabha events and muon 
pairs, the angle A0 in the x-y plane between the two secondary tracks was required to satisfy 
5° < < 175°. 

(6) To reduce the remaining background from two-photon events, each secondary track must have 
had a total momentum p > lOGeV/c. 

(7) To reduce the two-photon background, the invariant mass W of the two selected secondary 
tracks, which were assumed to be pions, was required to satisfy W > 5GeV/c^. 

For pair-produced smuons, selection efficiencies in the range of 60— 70 % were achieved. As selectrons 
would have been produced in the t-channel, leading to an enhanced production in the direction of 
the initial particle beams, selection efficiencies for selectrons range between 30 — 65 %. Efficiencies 
for chargino pair-production and decays to a smuon or selectron are in the range of 35-43 %. This 
efficiency is lower than for directly produced smuons or selectrons as additional visible energy is 
taken away by the neutrino from the chargino decay. In the whole data set no event was selected by 
the analysis, in agreement with the 0.13 it 0.03stat.i0.73syst. events expected from Standard Model 
sources. 

If the pair-produced NLSP is a stau or if the pair-produced chargino decays into a neutrino and 
the stau, the lepton in the event is a tau. As the tau decays to one or more charged particles plus a 
neutrino, additional energy is missing from the event. To take this into account, relaxed cuts were 
used for searches in the stau NLSP scenario. 

(6*) At least one of the secondary tracks must have had a momentum p > 5GeV/c. 

(7*) The invariant mass of the two secondary tracks was required to be at least 3GeV/c^. 

The efficiency for pair-produced staus ranges between 45 % and 50 %. For charginos in the stau 
NLSP scenario the efficiency ranges from 38 % up to 48 % and is a bit lower than for direct stau 
production because of additional missing energy from the neutrino from the chargino decay. After 
all cuts there were three events left in the data, while 3.79 it 0.57stat.il-14syst. events are expected 
from simulated Standard Model sources. 

3.2.2 Event selection for neutralino pair-production 

Pair-produced neutralinos are produced and decay via the processes e^e~ XiXi^ Xi ~^ ^R^ and 
in — > iG with the lepton a tau, muon or electron. Thus in this channel there are four leptons in 
the final state, two of which might have large impact parameters depending on the slepton NLSP 
lifetime. The cuts to search for such events were 

(1) events with at least two secondary tracks and at least two primary tracks were selected. No 
cut on the number of additional tracks was imposed. 

(2) For all secondary tracks the longitudinal impact parameter I^qI was required to satisfy l^^ol ^ 
40 cm to reduce cosmic particles, beam wall and beam gas events. 

(3) At least two of the primary tracks must have had a transverse momentum larger than 1 GeV/c. 
The transverse momenta of the primary tracks for the data and the simulated background 
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before this cut are shown in Figure 2. Apart from the poorly modeled very low-momentum 
region, mainly consisting of two-photon events, good agreement is found. After the cut, 158 
events remained in the data samples recorded at centre-of-mass energies between 189 GeV and 
209 GeV, while 186.4 events are expected from Standard Model sources. Additionally at least 
two primary tracks were required to have had hits in the silicon microvertex detector and the 
vertex chamber. 

(4) For Bhabha events the two primary tracks are usually back-to-back. To reject this background, 
the opening angle between all pairs of primary tracks ip in the x-y plane was required to satisfy 
if < 176°. 

(5) To separate the signal from the remaining two-photon background, the invariant mass W of 
all secondary tracks was required to fulfill W > b GeV/c^. In Figure 2 the distributions of the 
invariant mass before the cut are shown for data and simulated background. After the cut, 
12 events remained in the entire data sample, with 9.0 events expected from the background 
sources. 

(6) Secondary particles can also be produced in photon conversions. In this case their tracks often 
start in the middle of the jet chamber and have, due to the boost in the forward direction, 
relatively small impact parameters. For the signal, secondary particles arise from the decay 
of heavy SUSY particles, resulting in a much smaller boost. Thus to produce equally small 
impact parameters the flight lengths of the primary particles must be much shorter than for 
photons, and the secondary tracks start, in general, at radii smaller than those of the jet 
chamber. To veto such tracks from photon conversions the following cut was applied: for 
tracks with an impact parameter |(io| < 2 cm, the first hit wire in the jet chamber must have 
been measured at a radial distance to the primary interaction point of less than 40 cm. 

(7) Finally, to reduce remaining background from two-fermion and multi-hadronic events, the 
vector sum of the transverse momenta of all secondary tracks was required to exceed 3 GeV/c. 

Efficiencies of the selection are around 50 % in the slepton co-NLSP scenario and 35 % in 
the stau NLSP scenario. In total four events were left in the data in good agreement with 
5.71 ± 1.15stat.i0.29syst. events expected from simulated Standard Model background sources. 



3.2.3 Event selection for selectron or smuon pair-production (stau NLSP) 

The production and decay of selectrons or smuons following e'''e~ — > ^/j^^, (-R — > ^Xi; Xi ~^ '^^ 
f — > rG, leads to a final state with six leptons. In the event, two primary electrons or muons, two 
primary taus and finally two secondary taus and their decay products, which - depending on the 
NLSP lifetime - can have large impact parameters, are expected and the following cuts are used. 

(1) For a candidate event, at least two secondary and at least three primary tracks were required. 
No cut on the number of additional tracks was applied. 

(2) For all secondary tracks the longitudinal impact parameter \zq\ was required to satisfy |;Zo| !^ 
40 cm to reduce cosmic particles, beam wall and beam gas events. 

(3) At least two of the primary tracks must have had a transverse momentum larger than 1 GeV/c 
to reduce the background, mainly from two-photon events. Additionally at least two primary 
tracks must have had hits in the silicon microvertex detector and the vertex chamber. 
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(4) The invariant mass W calculated from all secondary tracks was required to be larger than 
5GeV/c^ to remove remaining background from two-photon events. 

(5) To veto photon conversions the radius of the innermost hit wire in the central jet chamber 
associated to tracks with {dol < 2cm was required to be smaller than 40cm, as described for 
cut (6) in Section 3.2.2. 

Selection efficiencies of 45% for both selectron and smuon pair-production are reached. In the entire 
data set, 13 events survived the selection, which is in good agreement with 13.72 it 1.74stat.±0.77syst. 
events that are expected from Standard Model background sources. 

3.2.4 Statistical and systematic uncertainties 

The statistical uncertainty on the number of expected background events due to the limited Monte 
Carlo statistics reach 50 % for single centre-of-mass energies mainly because of the limited statistics 
of the simulated two-photon events, which are expected to have high production cross-sections of 
the order of 10 nb. 

Most cuts applied in the search for tracks with large impact parameters rely on the tracking 
performance of the detector. If the resolution of the track parameters in the detector simulation is 
different from the one in the recorded data, this can lead to systematic effects on the number of 
expected background events. Therefore, to obtain a conservative estimation of the size of this effect, 
the resolution of the track parameters was degraded (smeared) by 10 % in the x-y plane and by 40 % 
in the z direction. With the smeared track parameters all analyses were repeated on the simulated 
background samples. The absolute differences to the background numbers without smearing were 
taken as the systematic uncertainties. For analyses with a very small number of expected SM events, 
the uncertainties can reach 100 % for single centre-of-mass energies. But in general systematic effects 
were found to be smaller than statistical uncertainties. 

3.3 Searches for charged NLSPs with medium hfetime (kinked tracks) 

The discovery signals for charged NLSPs with a lifetime in the range of about 10~9 s < r < 10"^ s 
are spectacular, with a kinked track in the tracking chambers. A kinked track is defined by the 
presence of a primary and at least one secondary track segment which intersect inside the sensitive 
volume of the central tracking chamber. In this analysis, all tracks reconstructed in the central 
tracking system are considered. 

Primary tracks: tracks originating at the interaction point of the event, traverse through 
the tracking detector volume with a minimum transverse momentum of llOMeV/c with no 
signal in the outer detectors. 

Secondary tracks: A minimum transverse momentum of llOMeV/c with |do| greater than 
2.5 cm and no associated hit in the silicon detector. 

Primary and secondary tracks can be combined into a kink candidate if there is an intersection 
point, the kink vertex, between the helices which describe the curvature of the primary and secondary 
tracks. No cut on the number of secondary tracks connected to a single primary track was made to 
keep sensitivity for heavy charged particles which possibly decay into a tau lepton. As a consequence 
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the search algorithms do not depend on the flavour of the slepton. Also, no selection cut was placed 
on the specific ionization energy loss dE/dx of the primary heavy charged particle, in order to keep 
the kinked track search independent from the searches for heavy stable charged particles. For the 
latter, the measurement of the dE/dx was the main tool. 

Depending on whether the NLSP was directly pair-produced or appeared in the decay chain of 
other SUSY particles, the kinked tracks would be accompanied by up to four additional charged 
leptons; thus, different search algorithms were designed and optimized for the different track multi- 
plicities in the event. 

3.3.1 Event selection for direct NLSP and chargino pair-production 

The processes e+e~ — > and e+e~ — > xfxi {xf — > Iri') followed by the delayed decay of the 

NLSP lead to a signature of two tracks in the event of which one or both can have a kink. To be 
selected as an event consistent with this topology the following criteria have to be satisfied: 

(1) At least one kink vertex in the event. 

(2) The sum of the numbers of primary plus secondary tracks must be less than or equal to 20 
and the number of good tracks should not exceed five. 

(3) There must not be more than two primary track candidates and the scalar sum of their 
transverse momenta pT had to be greater than 4GeV/c to reject events from two-photon 
processes. 

(4) The kink vertex had to be within the geometrical acceptance of the central jet chamber: 
r < 181.5cm and \z\ < (0.5 • [153.288cm +{r - 24.5cm) • tan(15°) • cos(7.5°)] - 1.0) . 

(5) The momentum of the secondary track had to satisfy p > IGeV/c and must not be parallel 
to a wire in the jet chamber to within 1° to remove combinatoric background and poorly 
measured tracks. 

(6) The transverse momentum of the primary track had to satisfy pt > I GeV/c. Figure 3 
(a) shows the distribution of the transverse momentum of the primary track. Unmodeled 
background events at momenta less than 1 GeV/c which are expected from two-photon events 
are removed. After the cut, 89 data events were left, consistent with 88.3 events expected from 
Standard Model sources. 

(7) The invariant mass Wqq of the primary particle was calculated using the primary and secondary 
track momenta and assuming the hypothesis of a massive particle decaying into two massless 
particles (one visible, and one invisible). Low-mass resonances were rejected by requiring this 
mass to be greater than 4GeV/c^. In Figure 3 (b) the distributions of the invariant mass for 
data and simulated background are shown. After the cut, 7 data events remained while 5.1 
events are expected from Standard Model sources. 

(8) Remaining background caused by shower electrons from Bhabha scattering which traverse the 
detector under small angles with respect to the beam axis were reduced by requiring that 
the angle between primary and secondary tracks be smaller than 172°. Also, a primary track 
having an angle smaller than 25° with respect to the beam axis must not point to a cluster in 
the electromagnetic calorimeter with an energy greater than 70 % of the beam energy. 
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(9) Hadronic interactions were suppressed by rejecting a kinked track with more than two sec- 
ondaries associated to its primary track. If there were two secondaries, the angle ^ between 
them was required to be smaller than 20° (cos C < 0.93) to keep sensitivity for decaying tau 
leptons. 

After all cuts 1.83ib0.42stat.i0.47syst. events were expected background for the entire data set. The 
background consists of events with hadronic interactions of particles with the detector material, 
mainly from two-photon processes. No data event passed all the cuts, which is consistent with the 
expectation from Standard Model sources. Typical efficiencies of the selection for pair-produced 
smuons and staus are around 45%, independent of the sparticle mass. For selectrons the efficiency 
decreases to 25 % at masses of 45 GeV/c^ as selectrons can be produced in the t-channel, which leads 
to a more forward production. 

3.3.2 Event selection for neutralino pair-production 

Searching for pair-produced neutralinos (e~'"e" — > XiXi? Xi ^ '^i''" or Xi ~^ ^R^) which decay into 
a lepton and the NLSP, the data were searched for topologies with a few tracks originating from 
the primary interaction point plus at least one kinked track. In the following list all cuts which are 
changed with respect to the search for pair-produced NLSPs (Section 3.3.1) are given. 

(2*) The sum of the numbers of primary plus secondary tracks must be less than or equal to 25 
and the number of good tracks should not exceed six. 

(3*) There must not be more than four primary track candidates and to reject events from two- 
photon processes the scalar sum of their has to be greater than 4GeV/c. 

(7*) The invariant mass Wqq of the primary particle has to be greater than 10 GeV/c^ to compensate 
for the relaxed cut on the number of tracks in the event. 

As for all other selections of events with kinked tracks the remaining expected background of 
1.05±0.35stat.±0.50 syst. events in this analysis consists of hadronic interactions from particles in 
two-photon events. One data event survived the complete selection which is consistent with the 
expectation. The efficiency for detecting neutralino pair-production followed by a delayed decay of 
the NLSP is typically around 40%, independent of the neutralino and NLSP mass. 

3.3.3 Event selection for selectron and smuon pair-production (stau NLSP) 

In the case of a cascade of pair-produced selectrons or smuons (e"'"e~ — > ^/j^/j, — > ixi: Xi ~^ 't^) 
to a stau with its decay inside the tracking chamber, the events contain a few tracks originating 
from the primary vertex plus at least one kink vertex. The cut on the number of tracks in the event 
has to be relaxed compared to the search for directly produced NLSPs. This is compensated by a 
tighter requirement on the invariant mass Wqq of the primary track. In the following list all cuts 
which are changed with respect to the search for direct NLSP production (Section 3.3.1) are given. 

(2**) The sum of the numbers of primary plus secondary tracks must be less than or equal to 25 
and the number of good tracks should not exceed ten. 
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(3**) There must not be more than eight primary track candidates, and to reject two-photon pro- 
cesses, the scalar sum of their transverse momenta pT should be greater than 4GeV/c. 

(7*) The invariant mass Wqq of the primary particle has to be greater than lOGeV/c^. 

After all cuts 2.76ib0.69stat.i0.51syst. events are expected from Standard Model sources, mainly 
hadronic interactions with the detector material. Three events survived the selection in the entire 
data set, consistent with the expectation. One surviving candidate is the same as the one selected 
by the search for pair-produced neutralinos (Section 3.3.2). The efficiency to detect the production 
and decay of pair-produced selectrons and smuons in the stau NLSP scenario is around 30 % for all 
slepton, neutralino and NLSP masses. 

3.3.4 Statistical and systematic uncertainties 

The expected Standard Model background for the kinked track search comes from hadronic interac- 
tions of particles with the detector material. All simulated processes with charged particles in the 
final state can contribute to it. Because of a high production cross-section, of the order of 10 nb, for 
two-photon events the number of simulated events of this type is not much higher than the number 
of data events expected. Thus selecting such an individual event contributes significantly to the 
number of expected background events and results in a significant statistical error, reaching up to 
100 % for a single centre-of-mass energy. 

The selections of events with kinked tracks are mainly based on tracking information. Thus 
systematic errors were estimated using the same method as the search for tracks with large impact 
parameters described in Section 3.2.4. The resolution of the parameters describing the tracks in the 
r-(j) plane in the simulated background and signal samples was degraded by 10 % and the resolution 
of the z parameters by 40%. The analyses were repeated on the smeared sample and differences 
with respect to the original results were counted as systematic errors. For a single centre-of-mass 
energy, differences up to 100 % were found. 

3.4 Searches for long-lived charged NLSPs 

If the slepton is the NLSP and the lifetime is such that the decay happens outside the detector 
volume the heavy charged particle has to traverse the central tracking detector. In this gaseous 
detector the measurement of the specific ionization energy loss dE/dx along the particle's track 
allows a powerful identification of signal events. 

3.4.1 Search for pair-produced stable charged NLSPs 

The search for pair-produced long-lived massive particles with a lifetime longer than 10~^ s is re- 
ported in [9]. The search is used here to select candidate events for pair-produced charged NLSPs, 
e''~e~ — > ^/j^^. The selection is primarily based on the precise measurement of the ionization energy 
loss {dE/dx) of charged particles in the jet chamber, using a data set corresponding to a total inte- 
grated luminosity of 693.1 pb~^ between centre-of-mass energies of 130 and 209 GeV. No candidate 
event was reported in the entire data set which is consistent with the expectation from Standard 
Model background of 1.1 events. 
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The typical efficiency of the search is about 95%; however, for particles with a (3j = p/m = 



Y 4^ — 1 between 0.75 and 1.70, the efficiency drops below 5%. Figure 4 shows that this effect is 
due to the fact that in this particular region, the dE/dx of tracks of heavy (SUSY) particles is 
very similar to that of light (Standard Model) particles. But as the LEP accelerator operated at 
different centre-of-mass energies, particles with a certain mass were in different efficiency regions for 
different centre-of-mass energies. 

3.4.2 Event selection for heavy stable charged particles in multi-track events 

If long-lived NLSPs are produced as secondary or tertiary particles from the decays of neutrali- 
nos, charginos or heavier sleptons, their tracks are no longer back-to-back and additional tracks or 
additional missing energy are expected. The analysis searching for directly pair-produced NLSPs 
(Section 3.4.1) cannot be applied and a different analysis was developed to search for tracks with 
anomalous ionization energy loss in events with more than two tracks or additional missing energy. 
For all signal topologies one rather general data selection is chosen where at least one track per 
event with an anomalous specific energy loss is required. 

The analysis was split into three steps: (1) preselection; (2) cuts on full event properties; (3) 
selection of high quality tracks with anomalous ionization energy loss. 

(1) Preselection 

(PI) Events with either more than 20 tracks or with a low total visible energy less than 0.10-^ 
were rejected to reduce the number of two-photon events. 

(P2) There must be at least one track in the event which has at least one hit in the precision vertex 
chamber, 20 hits used for dE/dx measurements and the specific energy loss should be higher 
or lower than that expected from SM particles {dE/dx > llkeV/cm or dE/dx < 9keV/cm 
for a track momentum p > lOGeV/c). 

(P3) Non-simulated backgrounds like cosmic muons and beam gas interactions are removed by 
requiring that the distance between the beam axis and the track at the point of closest approach 
had to be smaller than 1.5 cm in the x-y plane and smaller than 10 cm along the z coordinate. 
If there were hits in the time-of-flight barrel scintillators the measured time of the closest 
barrel hit was allowed to differ by at most 10 ns from the time expected from a physics event. 

(2) Event Analysis 

(1) To reduce two-photon background events, which deposit a large fraction of energy at small an- 
gles with respect to the beam axis, the maximum energy deposited in the forward calorimeters 
must be less than 5 GeV. 

(2) The total measured relativistic invariant mass Mvis (visible mass) must satisfy Mvis > 
lOGeV/c^ to reduce further the number of two-photon events. 

(3) The visible energy was required to satisfy 0.15 < E'vis/v^ < 1.10, to reject two-photon and 
two-fermion events. 
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(4) The maximum electromagnetic energy -Eecal was required to satisfy E-^q^i^/ ^/s < 0.3 to 
separate the signal from two-fermion processes. 

(5) If the sum over all track momenta in an event P = |P| = |X]Ptrack| was greater than 10 GeV/c, 
the z component of the sum Pz must satisfy \Pz/P\ < 0.9. This cut reduced mainly Bhabha 
scattering events with one electron escaping through the beam pipe carrying away a significant 
fraction of the total momentum of the event. 

In the data set recorded between 189 GeV and 209 GeV, 572 data events survived the requirements 
on the event properties while 537.04 events are expected from Standard Model sources. Figure 4 (a) 
shows the distribution of the variable \Pz/P\ used in cut (5) together with the distribution of an 
expected signal for pair-produced staus. 

(3) Anomalous dE/dx selection 

(6) In case the track with momentum p is associated to a cluster in the electromagnetic calorimeter 
with an energy Ecai, the condition Ec^i/p < 0.15 has to be satisfied. This cut further reduces 
Bhabha scattering events, four-fermion and two-photon background. 

(7) To avoid events with converted photons - where the electron and positron tracks are unresolved 
by the tracking detectors, hence give a single track with a high dE/dx - the track must have 
at least one hit in the high precision vertex chamber. 

(8) The track had to satisfy tightened quality criteria. At least 40 hits should be used for dE/dx 
measurement and the error on the momentum must be smaller than lOGeV/c. 

(9) The track with momentum p must feature a high or low specific ionization energy loss dE/dx: 

if p > a • {dE/dx - b) then dE/dx > 12.0keV/cm (a = ^ cm, b = 181.5 keV/cm) or 
if p > 52GeV/c then dE/dx < 8.2keV/cm 
The cut separated signal topologies from Standard Model background. 

(10) Finally, to separate the signal from Standard Model events with taus, no other track was 
allowed in a cone of 20° around the track with the anomalous dE/dx. 

Figure 4 (b) shows the distribution of the candidate events before cut (9) as a function of the 
track momentum and the specific ionization energy loss. The regions in which candidate events are 
expected are indicated. Out of 448 candidate events (419.33 events expected from Standard Model 
sources), no data event survived selection cut (9) while 0.89 events are expected. 

The typical efficiency of the search ranges between 80 % and 95 % for masses of the heavy 
particle near the beam energy and between 20% and 40% for particles with a (3j larger than 1.5. 
The efficiency drops below 5 % only in a small region of (3^, where the ionization energy loss for 
massive particles is very similar to the energy loss of Standard Model particles, as shown in Figure 4. 
In the data, recorded at = 130 GeV~209 GeV and corresponding to an integrated luminosity of 
C = 632.1 pb~^, no event survived after all cuts. This is compatible with 0.78±0.38stat.±0.10syst. 
events which are expected from the Standard Model sources. 
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3.4.2.1 Systematic uncertainties 

The main tool used in this analysis is the information from particle tracks. Therefore, the 
main systematic errors of the event selection arise from uncertainties in the track quality and the 
modeling of the specific energy loss in the simulation. 

The systematic error arising from the track measurement was evaluated by smearing simulated 
track resolution by 5 % for r-(j) and 20 % for z. Repeating the entire analysis with modified track 
parameters for signal and background determines the first contribution to the systematic error 
arising from the central tracking system. The second contribution to the systematic error is due 
to uncertainties in the modeling of the measured dE/dx value. These uncertainties were evaluated 
by comparing muon pairs, Bhabha electrons, charged pions, and kaons in simulation and data, and 
found to be 10 % of the measurement error (J^e/Ax- The contribution to the overall systematic error 
was determined by redoing the analysis of signal and background with the dE/dx value of each 
track replaced by a value {dE/dx ± (t^e/Ax)- 

The studies show that the number of expected background events may vary by 10 % from the 
track smearing and by 5.0% due to the dE/dx uncertainties. 

4 Searches in the neutrahno NLSP scenario 

In this section the analyses to search for indirect NLSP production in the neutralino NLSP scenario 
are described. To be sensitive to all neutralino lifetimes, searches for the production of sleptons and 
charginos were performed. Section 4.1 describes the analyses selecting event signatures expected 
for short-lived neutralinos with lifetimes up to 10~®s. In Section 4.2 searches for neutralinos with a 
long lifetime are presented. A summary of the results for searches not described in other papers is 
listed in Table 3. 

4.1 Searches for short-lived neutrahnos 

All signatures of SUSY particles in the scenario with the neutralino being the NLSP have in common 
that, for short and medium lifetimes up to 10~® s, photons from the neutralino decay would be 
measured by the detector. For longer NLSP lifetimes, the decay happens outside the detector and 
typical signatures contain a significant amount of missing energy. 

4.1.1 Search for events with photons and missing energy 

The search for pair-produced neutralinos, decaying in the detector volume, selects events with two 
photons and large missing energy. The main background for the search is the Standard Model 
process e"'"e~ uv + photons. The selection and the result of the search are described in [10]. No 
evidence is observed for new physics contributing to the expected final state. 

4.1.2 Selection of events with isolated photons, missing energy, leptons and jets 

If the lightest neutralino is the short-lived NLSP, the production and decay of scalar leptons (e"'"e~ — > 
^~R^R^ ~^ ^Xi) or charginos (e+e" — > xtxii xt ~^ W^Xi) lead to final states with missing 
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energy, up to two photons plus leptons and/or jets (7(7)-^ plus qiqjqkqi) QiClj^^ or For 
each final state, the photon energy and the missing energy depend strongly on the mass difference 
AM (=Mnnlsp — -^nlsp) between the pair-produced particle (the NNLSP) and the neutralino (the 
NLSP). Thus a set of analyses was designed, using a common preselection and followed by specific 
series of cuts optimized for the different AM regions: 

(A) photons plus leptons and missing energy; 

(B) photons plus jets and missing energy 
(BO) small AM (< 20 GeV/c^), 

(Bl) medium AM (20 < AM < Mnnlsp - 30 GeV/c^), 
(B2) large AM (> Mnnlsp - 30 GeV/c^); 

(C) photons plus jets plus a lepton and missing energy 

(CO) small AM (< 20 GeV/c^) where analysis (BO) is used, 
(CI) medium AM (20 < AM < Mnnlsp - 30 GeV/c^), 
(C2) large AM(> Mnnlsp - 30 GeV/c^). 

The first part of the preselection required well-contained events from e~''e~ collisions and is identical 
to that used in the long-lived neutralino NLSP chargino search given in Section 4.2.2 and documented 
in [11]. The lepton and jet identification used here is also identical to that in [11]. All events with 
more than six good tracks from the interaction point and not originating from photon conversions 
were classified as candidates for selection (B) or (C), events with fewer tracks as candidates for 
selection (A). 

The second part of the preselection picks up events with isolated photons and tightens the 
topological requirements. 

(PI) There have to be at least two isolated photons with: 

- an electromagnetic cluster with an energy greater than 3 GeV; 

- the angle of the cluster with respect to the beam satisfying | cos9\ < 0.9; 

- in a 15° half-cone centred on the cluster the sum of track momenta and additional 
electromagnetic energies not exceeding 2 GeV. 

(P2) The angle ^miss of the missing momentum had to satisfy | cos^missj < 0.95. 

(P3) The transverse momentum of the event had to be greater than 3 % of the centre-of-mass energy. 

(P4) The visible energy in the event had to satisfy 

- 0.2 • ^ < E^is/^/s < 0.9 • ^ for analysis A, 

- 0.25 • ^/s < E^is/V^ < 0-85 • ^/s for analysis BO and 

- 0.4 • ^/s < E^is/^/s < 0.95 • ^/s for the other analyses. 

In Figure 5 (a) the distribution of the visible energy in the event is shown for the data recorded 
at a centre-of-mass energy of 206 GeV together with the expected background. After cut (P4) 
of analysis A, 354 data events remained in the sample with 368.48 events expected in the 
background. Plot (b) shows the same distribution for a possible signal of stau pair-production 
at three different mass combinations of the stau and the neutralino. 
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(P5) The acoplanarity angle (pacop, which is defined as the supplementary angle in a plane perpen- 
dicular to the beam direction between the two vectors which are obtained by summing up 
particle momenta after splitting the event into two jets using the Durham algorithm, has to 
satisfy (/)acop > 5°. 

After the preselection, the cuts were optimized depending on the mass difference between the NNLSP 
and the NLSP as well as the number of leptons expected in the event. 

Selection (A): photons plus two leptons and missing energy 

Events with two photons, two leptons and missing energy are expected from the production 
and decay of sleptons (e+e~ — > i^i]^, — > ^x?, Xi ~^ 7G) and charginos (e"''e~ XiXii 
Xi — > W^Xij Xi 7G) with both W bosons decaying leptonically. The main expected back- 
grounds from Standard Model sources are radiative return events to the Z boson, with TP — > 
accompanied by initial or final state radiation. To suppress this type of background, cuts on the 
photon energies are sufficient. The analysis does not distinguish between electrons, muons or taus 
since for small mass differences AM the lepton energy is small and they cannot be identified 
efficiently. 

(A-1) The energy of the most energetic photon had to be between 10 GeV and 90 GeV, and 
the energy of the second photon greater than 5 GeV. Almost all radiative return events were 
rejected by the requirement on the second most energetic photon. In Figure 5 (c) the energy 
distribution of the most energetic photon is shown for the data recorded at a centre-of-mass 
energy of 206 GeV together with the expected background. After the cut, 74 data events 
remained in the sample with 70.63 events expected in the background. Plot (d) shows the same 
distribution for a possible signal of stau pair-production at three different mass combinations 
of the stau and the neutralino. 

(A-2) Surviving events with final state radiation often have the photon near a particle track. To 
reject them the transverse momentum of the photon with respect to the nearest track had to 
be larger than lOGeV/c (5GeV/c) for the first (second) photon. 

In data recorded at centre-of-mass energies between 189 GeV and 209 GeV, 3 events survived the 
selection with 7.74ib0.47stat.il-63syst. events expected from Standard Model sources. 

Selection (BO): photons plus jets and missing energy for small AM 

In the case of production and decay of charginos (e+e~ — > Xi^xr' ~^ W^Xi; Xi ~^ 7*5) 
with both W bosons decaying into two jets for small mass differences between the chargino and the 
neutralino NLSP, the jets contain only a small amount of energy. The jets might not be measured 
well, thus the selection depends on the two most energetic photons only. 

(BO-1) The energy of the most energetic photon has to be greater than 20 GeV and the energy of the 
second photon greater than 10 GeV. 

(BO-2) The visible energy in the event, excluding the two most energetic photons, has to be less than 
50 GeV. 

After the selection there is no event left in the entire data set with 0.69ib0.09stat.i0.14syst. events 
expected from Standard Model sources. 
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Selection (Bl): photons plus jets and missing energy for medium AM 

In case of a higher mass difference between the chargino and the neutralino, the four jets of 
the hadronically decaying W bosons have significant energy, and there is no prompt lepton in the 
events. The jets exclude the identified photons. 

(Bl-1) The number of identified, isolated leptons had to be zero. The lepton was defined to be isolated 
if there was not more than 2 GeV additional energy in a cone with a half-angle of 15° around 
the lepton. 

(Bl-2) To reject events with final state radiation, the transverse momentum of the photon with respect 
to the nearest track had to be larger than lOGeV/c (5GeV/c) for the first (second) photon. 

Different additional requirements have to be fulfilled, depending on the number of jets in the event. 

(Bl-3) 4-jets: the angle between the two reconstructed W bosons had to be less than 150°, with the 
two bosons (jet pairs) being selected by minimizing the difference between the jet masses and 
the W boson mass ((M^ - Mwf + {M^i - Mwf)- 

3-jets: the two jets with an invariant mass closest to the mass of the W boson were paired 
together. The angle between the W and the remaining jet had to be smaller than 150°. 

2-jets: the sum of the opening angles between each jet and the most energetic photon had to 
be less than 356°. 

After the selection 6 data events survived the cuts of the medium AM analysis (Bl) with 
8.59ib0.86stat.=tl-80syst. events expected from Standard Model sources. 

Selection (B2): photons plus jets and missing energy for large AM 

Because of the large mass difference between the chargino and the neutralino, expected sig- 
nal events have less energetic photons and only a small missing momentum. After the selection, 
the remaining backgrounds are mainly from the processes e"'"e~ W'^W" and e"'"e~ — > TP^ with 
additional photons. 

(B2-1) The number of identified, isolated leptons had to be zero. 

(B2-2) To reject events with final state radiation, the transverse momentum of the photon with respect 
to the nearest track had to be larger than 5GeV/c (3GeV/c) for the first (second) photon. 

(B2-3) The event was forced into four jets and the invariant mass cut of the Durham jet finding 
algorithm kx had to satisfy kx > 4GeV/c^. 

After the selection, 3 data events survived with 2.67ib0.26stat. events expected from Standard Model 
processes. 



25 



Selection (CI): photons plus jets plus one lepton and missing energy for medium AM 

If one W boson decays into a lepton and a neutrino, the other hadronically, two photons, 
two jets and one lepton with lower energies are expected in the event. The lepton can be identified 
and candidate events had to satisfy the following constraints. 

(Cl-1) There had to be at least one isolated lepton in the event. 

(Cl-2) To reject events with final state radiation, the transverse momentum of the photon with respect 
to the nearest track had to be larger than lOGeV/c (5GeV/c) for the first (second) photon. 

(Cl-3) The lepton energy had to be less than 40 GeV and the invariant mass of the two jets less than 
70 GeV/c^. This reduced the background from W"^W~ — > qqiv where jets and lepton are 
more energetic compared to expected signal events. 

One event survived the selection cuts for the medium AM analysis, with 0.62ib0.11stat. iO.lSsyst. 
events expected from Standard Model sources. 

Selection (C2): photons plus jets plus one lepton and missing energy for large AM 

In case of a large mass difference between the chargino and the neutralino the jets and lep- 
ton have more energy. The main expected background results from the decay of pair-produced W 
bosons, W~'^W~ — > qqii^i, with additional photon(s) from the initial and(or) final state radiation. 
The signal topology is similar to the Standard Model process, but it can be separated by requiring 
two isolated photons and unbalanced W bosons. 

(C2-1) There had to be at least one isolated lepton in the event. 

(C2-2) The transverse momentum of the photon with respect to the nearest track or jet had to be 
larger than 5GeV/c (3GeV/c) for the first (second) photon. 

(C2-3) The angle between the two W bosons, where one W is reconstructed from the two jets and 
the other by the most energetic lepton and missing momentum, had to be less than 170°. 

In the data sample, two events survived the selection cuts, with 2.94ib0.31stat. events expected from 
Standard Model processes. 

4.1.2.1 Systematic uncertainties 

The uncertainty associated with the photon isolation requirement is the dominant system- 
atic error in these selections. It originates from the modeling of parton-level photon emission, jet 
fragmentation and detector simulation. Since it is difficult to estimate an error on each source, a 
total error is obtained inclusively by comparing data with Monte Carlo simulation using several 
cross-check samples. The studies are limited by statistics and a systematic error, which is at most 
20%, is taken into account for all selections. 

Uncertainty on the modeling of the other cut variables for simulated background and signal 
was also investigated. This was done by shifting the cut value within the possible error. This is 
negligible for signal events but not for the background, especially at large AM, because the topology 
is similar to Standard Model events and some cuts are set near to the peak of their distributions. 
Uncertainties from 3 % for analyses (A) and (BO) up to 6 % for the other analyses are reached. 
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4.2 Searches for long-lived neutralinos 



In the case that the neutrahno is the NLSP and has a hfetime longer than 10~ s, it usually decays 
outside the detector volume. The topologies are similar to searches for supersymmetric particles in 
gravity-mediated SUSY breaking scenarios in which the stable neutralino is the lightest supersym- 
metric particle. 

4.2.1 Search for di-lepton events with missing transverse momentum 

If the pair-produced sleptons decay to a lepton and a long-lived neutralino, the signature of possible 
candidate events is two leptons with missing transverse momentum. The search for this type of 
events is reported in [7] and a short description can be found in Section 3.1.1. 

The search for di-lepton events with missing transverse momentum is also sensitive to the pair- 
production of charginos. Here the chargino can decay either to a slepton and a neutrino, followed by 
the prompt decay of the slepton to a lepton and a long-lived neutralino, or the chargino can decay 
to a neutralino and a W boson with both W bosons decaying leptonically. 

4.2.2 Search for charginos decaying into a neutralino and a W boson 

Approximately 438 pb^^ of data recorded at centre-of-mass energies of 192-209 GeV were analyzed 
to search for evidence of chargino pair-production. No evidence for a signal was observed and the 
results of the search are reported in [11]. Search channels looking for charginos decaying into a 
neutralino and a W boson were used for the GMSB searches described in this paper. All decays 
of the W boson, the hadronic or leptonic decay, leading to topologies containing jets and missing 
energy, or jets with a lepton and missing energy were considered. 

In case both W bosons decay leptonically, the results of the searches reported in [7], Section 4.2.1 
were applied. 

5 Analysis combination strategy 

For all sparticle production and decay channels the results from several analyses at several centre- 
of-mass energies have to be combined. Although the different analyses are sensitive to distinct 
topologies, some correlation, an overlap between selected data, background and signal events, is 
expected. To treat the overlaps between the various analyses properly, all used the same signal 
Monte Carlo samples, which were generated for each search channel, as well as the same samples of 
the simulated Standard Model background. This allows, via an event-by-event comparison of the 
selected signal, expected background and selected data events, the determination of the overlaps 
between all combinations of the analyses used in each channel. 

The selection efficiencies were split into exclusive selection efficiencies and overlap efficiencies: 
the efficiency for events selected by only one analysis are denoted as exclusive selection efficiencies; 
events that are selected by two or more analyses are described by overlap efficiencies. For example, 
for the case of three analyses, there are, in total, three exclusive selection efficiencies and four overlap 
efficiencies (three for the overlaps between each two analyses and one for the overlap among all three 
analyses). 
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Sparticle 
production 
(NLSP) 


NLSP 
lifetime 


< s/s> = 
188.7 GeV 
(168 pb-i) 


< \/s > = 
191.6 GeV 
(30pb-i) 


< \/s > = 
195.5 GeV 
(76 pb-i) 


< ^ >= 
199.5 GeV 

(78pb-i) 


<^>= 
201.6 GeV 
(38pb-i) 


< \/s >= 
205.1 GeV 
(74pb-i) 


< \/i >= 
206.7 GeV 
(120pb-i) 


< \/s >= 
208.1 GeV 
(8pb-i) 


Reference 


(h) 


(1) short 


l/0.93±0.56 


0/0.11±0.06 


1/0.25±0.13 


0/0.49±0.36 


0/0.18±0.28 


0/0.56±0.24 


1/1.13±1.02 


0/0.14±0.04 


3.2.1 


(2) medium 


0/0.71±0.31 


0/0.09±0.12 


0/0.22±0.33 


0/0.29±0.28 


0/0.17±0.18 


0/0.13±0.14 


0/0.21±0.22 


0/0.01±0.02 


3.3.1 


xlxl 
in) 


(3) zero 


1/1.06±0.45 


0/0.15±0.11 


0/0.38±0.18 


0/0.31±0.17 


0/0.30±0.15 


0/0.24±0.15 


l/0.38±0.22 


0/0.03±0.04 


3.1.2 


(4) short 


1/1.20±0.33 


l/0.23±0.07 


0/0.57±0.31 


0/1.25±0.90 


l/0.23±0.05 


0/0.91±0.31 


1/1.24±0.53 


0/0.08±0.04 


3.2.2 


(5) medium 


0/0.69±0.52 


0/0.00±0.01 


0/0.01±0.01 


l/0.06±0.04 


0/0.04±0.03 


0/0.09±0.18 


0/0.15±0.28 


0/0.01±0.02 


3.3.2 


(6) long 


0/0.07±0.05 




0/0.08±0.02 


0/0.33±0.22 


0/0.04±0.01 


0/0.06±0.02 


0/0.11±0.04 


0/0.01±0.00 


3.4.2 


xtxi 
(n) 


short NLSP lifetime as for (1), medium It. as for (2), long It. as for (6) 


in) 


(7) zero 


2/1.59±0.60 


0/0.22±0.13 


0/0.56±0.23 


l/0.46±0.20 


0/0.46±0.19 


l/0.35±0.20 


1/0.57±0.31 


0/0.05±0.05 


3.1.2 


(8) short 


4/3.68±1.08 


1/0.64±0.19 


1/1.00±0.37 


2/1.91±0.92 


2/0.62±0.16 


l/2.23±0.60 


2/3.42±1.01 


0/0.22±0.066 


3.2.3 


(9) medium 


1/1.34±0.67 


0/0.13±0.09 


0/0.33±0.23 


1/0.41±0.30 


0/0.25±0.13 


1/0.11±0.18 


0/0.18±0.28 


0/0.01±0.02 


3.3.3 


long NLSP lifetime as for (6) 


in) 


zero NLSP lifetime as for (7), short It. as for (8), medium It. as for (9), long It. as for (6) 


*^R.^R. 

(Cb) 


(10) short 


0/0.05±0.03 


0/0.01±0.01 


0/0.01±0.01 


0/0.02±0.01 


0/0.01±0.01 


0/0.01±0.01 


0/0.01±0.73 


0/0.01±0.01 


3.2.1 


medium NLSP lifetime as for (2) 


xlxl 
{iR) 


zero NLSP lifetime as for (3), short It. as for (4), medium It. as for (5), long It. as for (6) 


xtxi 


short NLSP lifetime as for (10), medium It. as for (2), long It. as for (6) 


^R^R 

[(11) only] 
and 

xtxi 

(X?) 
[(11) + (12)] 


(11 A) short 


0/1.97±0.45 


0/0.43±0.10 


0/1.02±0.24 


0/1.10±0.25 


0/0.45±0.11 


2/1.03±0.31 


1/1.62±0.48 


0/0.12±0.04 


4.1.2 


(12 BO) short 
(12 Bl) short 
(12 B2) short 
(12 CI) short 
(12 C2) short 


0/0.13±0.05 
0/1.71±0.40 
0/0.60±0.15 
0/0.12±0.08 
2/0.60±0.15 


0/0.03±0.01 
0/0.45±0.12 
0/0.12±0.04 
0/0.03±0.01 
0/0.09±0.03 


0/0.03±0.01 
l/0.99±0.26 
0/0.27±0.09 
0/0.06±0.03 
0/0.25±0.07 


0/0.04±0.01 
3/1.17±0.29 
l/0.30±0.09 
0/0.04±0.01 
0/0.22±0.06 


0/0.03±0.01 
0/0.43±0.18 
1/0.12±0.04 
0/0.05±0.02 
0/0.16±0.04 


0/0.16±0.05 
0/1.43±0.48 
0/0.47±0.16 
0/0.12±0.04 
0/0.60±0.20 


0/0.25±0.09 
2/2.25±0.75 
l/0.74±0.25 
1/0.19±0.06 
0/0.95±0.32 


0/0.02±0.00 
0/0.16±0.41 
0/0.05±0.01 
0/0.01±0.00 
0/0.07±0.02 


4.1.2 
4.1.2 
4.1.2 
4.1.2 
4.1.2 



Table 3: Number of observed and expected events for all searches described in this paper. The error on the expected number of events 
include both statistical and systematic effects. In the interpretation, results of analysis (6) at lower centre-of mass energies are used 
in addition. Here the numbers of observed and expected events are 0/0.01±0.00 at ^/s = 130 GeV, 0/0.01±0.00 at ^/s = 161 GeV, 
0/O.OlibO.OO at ^/s = 171 GeV and 0/0.05ib0.01 at ^/s = 183 GeV. The numbers given in brackets in the first row of the table give 
the average recorded luminosity in each energy bin. 



For all search channels the signal efficiencies, i.e. exclusive selection efficiencies and significant 
overlap efficiencies, were determined at each generated mass-lifetime point of the signal Monte Carlo 
samples. For data as well as for the expected background events the overlaps between the different 
analyses were investigated. The overlap for the background was found to be negligible (< 0.01 
events) in all slepton NLSP channels. No data event was selected by more than one analysis. The 
various overlap selections were treated as separate search channels when combining them with the 
exclusive selections in the limit calculations. 

In order to avoid generating and processing excessive amounts of signal Monte Carlo samples 
and to achieve a good description of the efficiencies over the whole mass and lifetime range an 
interpolating function was determined to calculate the efficiencies at any given mass, lifetime point 
and centre-of-mass energy. This is particularly important in the case that the NLSP occurs at the 
end of a cascade decay. In this case, it is impossible to simulate all possible mass combinations of 
the SUSY particles, hence the efficiency of the analyses has to be determined by the kinematics of 
the NLSP. In the next section such a fit function is discussed. 



5.1 The efficiency function 

The general ansatz for the interpolating function uses the following assumptions: 

• The decay of the NLSP is described by exp(— where t is the time since the NLSP production 
and r is the NLSP lifetime. 

• The detector to measure the decay is a sphere with a fiducial volume ranging from a radius 
/start to a radius /stop- 

• The detection efficiency of the detector ea and the efficiency of the analysis e^, (overall efficiency 
Co = Crf ® Ca) are constant within the bounds /start and /stop- 

• Initial and final state radiation are neglected, which is a good assumption for heavy particles 
near the kinematic limit. 



The probability of detecting a particle is a function of its mean decay length /meam which is a 
function of its lifetime, r, and its boost, Pj. For pair production, (3^ is a function of the particle mass, 
TTT-i, and centre-of-mass energy ^/s. Then the probability of an analysis with an overall efficiency eo 
detecting a particle decay inside a detector with the radii /start and /stop is given by 

/start \ / /stop \ /o\ 

~'i — ]""''Pr/ — j ' 

'mean / \ "rncan / . 

with /mean, the decay length of the particle, defined as 



/mean = Cr/3i7i = cry'^^ - 1. (3) 

The values Csingie; 4tart and /stop are effective values which are determined by a fit of the efficiency 
function (2) to the efficiencies obtained for the signal Monte Carlo samples, treating each channel 
separately. Depending on whether an analysis accepts events where one or both pair-produced 
particles have to decay in the detector, the event selection efficiency eevent is given by 

eevent = 1 - [1 - Csingic]^ or eevent = Csingio respectively. 



^single ^0 



exp 
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These formulae are used to describe the selection efficiency of the search analysis as a function of 
the NLSP (3^ over a wide range of masses, lifetimes and production energies with three parameters, 
^0) 4tart and /stop- The input variables are the mass of the decaying particle, its lifetime and the 
centre-of-mass energy. For the calculation of the efficiency of the zero lifetime search, the parameter 
/start can be set to zero as the sensitive volume starts at the interaction point. For the heavy charged 
stable particle search, the parameter /stop can be set to infinity. 

An example of the efficiencies described by the function is given in Figure 6 for stau pair- 
production at ^/s = 208 GeV. Here the search for acoplanar leptons, for tracks with large impact 
parameters or kinks, the search for particles with anomalous specific ionization energy loss plus all 
overlap analyses are included. 

5.1.1 The efficiency function for NLSP in cascade decays 

If the particle Ti with lifetime r is not pair-produced but a decay product of another pair-produced 
particle Tq with mass mo and energy Eo = ^, the average decay length /mean has to be replaced by 

/mean = l^lllCT = — CT . 

mi 

Here \pi \ , the average momentum of Ti , is calculated from the two-body decay of To into Ti and a 
second decay product. 

5.1.2 Selection efficiency for long-lived, charged NLSP from cascade decays 

The use of the average monientum in the cajlculation of the NLSP decciy length /mean 

as described 

in the previous section is only justified if the efficiency, as a function of the NLSP /37, is a smooth 
function with moderate variation of the order of a few percent. As the efficiency of searches for 
long-lived, heavy particles from cascade decays drops drastically from about 80 % to below 5 % for 
a {3^ of the particle between 0.7 and 1.7, a different method has to be applied. Instead of using 
the average NLSP momentum \p\/m, the kinematically possible (3^ spectrum after the decay of a 
primary particle is used, following the ansatz 

e = y eo(/3i7i,^272) • 5(/5i7i> /3272) d/3i7id/3272 • 

Here, eo is the efficiency to select a signal event expressed as a function of the j3^ values of the two 
possible NLSPs. It is determined using simulated signal events. The normalized (3^ spectrum of 
the heavy stable charged particles is given by the function g. The spectrum is calculated assuming 
a flat energy spectrum of the stable particles in the kinematically accessible range and neglecting 
ISR. 

5.1.3 Errors on the interpolated efficiencies 

In addition to the systematic uncertainties on the selections described in the analysis sections, an 
additional source of the systematic uncertainty for the signal selection results from the interpola- 
tion of the efficiencies to different mass and lifetimes. Its size was estimated by recalculating the 
parameters of the efficiency function dropping the information of one simulated signal sample and 
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comparing the interpolated efficiency with the achieved efficiency at this point. Deviations were 
found to be less than 3 %. Finally the efficiency function was varied by the errors on the fitted 
parameters eo, ^start and /stop- Here an uncertainty of less than 1 % was found. The total systematic 
error assigned to the efficiency function is calculated by adding the three terms in quadrature. 

The efficiency function for long-lived, charged NLSPs in cascade decays, Section 5.1.2, has regions 
of rapid variation and an additional systematic error is assigned. In the region of NLSP [3^ between 
0.75 and 1.70 the efficiency is varying rapidly to small values, and relative uncertainties up to 43 % on 
the modelling and 31 % on the interpolation of the function were found. In the high efficiency region, 
the two contributions were determined to be 11% on the modelling and 8% on the interpolation. 
These dominant effects were included in the calculation of the limits on the production cross-section 
of SUSY particles. 

6 Constraints on particle production cross-sections 

The absence of any significant excess of events in the data compared to the expected number of 
Standard Model events in any of the search analyses can be translated into limits on the production 
cross-section of the SUSY particles at the 95 % confidence level. The limits were calculated using 
the program described in [31], which incorporates statistical and systematic uncertainties on the 
efficiency and the expected background as well as the uncertainties on the luminosity of the analyzed 
data samples and uncertainties on the expected cross-sections into the limits using a numerical 
convolution technique. The results for each analysis were given at eight different bins of the centre- 
of-mass energy (see Table 3). For all analyses, systematic and statistical uncertainties on the signal 
efficiency and on the background expectation as well as uncertainties on the recorded luminosity 
were included. 

Deriving cross-section limits at a single centre-of-mass energy requires knowledge of the cross- 
section evolution with ^/s and, in general, this evolution depends on the details of the SUSY model. 
For direct NLSP pair-production, the GMSB model database described in Section 7 is scanned, and 
the cross-section evolution for each point in the database is used to calculate the cross-section limit 
at y/s = 208 GeV for each mass and lifetime. The maximum excluded cross-section limit is chosen 
as the "model independent" limit for this mass and lifetime. For "lifetime independent" exclusions, 
the maximum cross-section limit valid for all lifetimes is chosen. 

In the channels with cascade decays to the NLSP, this procedure does not work because there 
are points in the GMSB model with vanishing NNLSP production cross-sections and therefore the 
evolution is not defined. Instead, in these cascade channels, for spin 1/2 supersymmetric particles 
a, (5/s dependence of the cross-section is assumed, while for scalars a dependence is used. 

The cross-section limits quoted in the cascade channels were evaluated for the "worst-case" of all 
intermediate particle masses, as well as NLSP lifetime. 

6.1 The stau or slepton as the NLSP 

Here the results of searches for signatures expected if the NLSP is the lightest slepton are described. 
The sleptons are either produced directly or appear as decay products of heavier pair-produced 
SUSY particles. 
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6.1.1 Cross-section limits for direct NLSP (stau or slepton) pair-production 

To obtain constraints on the production cross-sections for slepton pair-production in the slepton 
co-NLSP scenario and stau pair-production in the stau NLSP scenario the search for acoplanar 
leptons (Section 3.1.1), the search for tracks with large impact parameters (Section 3.2.1), the search 
for kinked tracks (Section 3.3.1) and the search for heavy stable charged particles (Section 3.4.1) 
were combined. A significant overlap of up to 40 % in the signal efficiency was found between the 
acoplanar lepton and large impact parameter searches. The overlap between the acoplanar lepton 
and kink searches is about 4%, between the large impact parameter and kink searches about 7%, 
and between the kink and heavy stable charged particle searches up to 20%. Examples of the 
exclusive and overlap signal efficiencies can be seen in Fig. 6 for the stau searches. There is no data 
event selected in common, while for the expected background an overlap of 0.01 events was found 
between the large impact parameter and acoplanar lepton searches in the selectron and the stau 
channels. 

The resulting limits at 95% C.L. at ^/s = 208 GeV are given in Figure 7 for staus in the 
stau NLSP scenario (a) and for staus (b), smuons (c) and selectrons (d) in the slepton co-NLSP 
scenario. For all ffavours the weakest exclusion is found in the very short lifetime region (r = 
10~^^ — lO"^'^ s), due to the irreducible background from W~''W~ production in the acoplanar lepton 
topology, whereas all other regions are almost free of expected background. Also, the smaller 
sensitivity for masses around 60GeV/c^ at long lifetimes, due to the loss of sensitivity at d£^/dx 
band crossings as described in Section 3.4.1, is visible in all channels. Because of the highest 
selection efficiencies, the best constraints are obtained for smuons with an upper limit of 0.05 pb 
in most of the plane, apart from the region with a smuon mass above lOOGeV/c^, close to the 
kinematic limit. For selectrons the constraints are slightly weaker. Here cross-sections larger than 
0.1 pb can be excluded in most of the selectron mass-lifetime plane. The difference between smuons 
and selectrons is explained by the fact that selectrons can be produced in the t-channel, which leads 
to a more forward production and thus to lower efficiencies. For stau pair-production, cross-sections 
larger than 0.1 pb are excluded. The stau limits in the slepton co-NLSP scenario and in the stau 
NLSP scenario are rather similar, with a small difference due to the slightly different theoretical 
cross-sections that are used for the combination of the results at various centre-of-mass energies. 

6.1.2 Cross-section limits for neutralino pair-production 

Two scenarios exist in the theory: in the stau NLSP scenario the neutralino decays with 100 % 
branching ratio to fir, while in the slepton co-NLSP scenario it decays with equal branching fractions 
to all flavours. Both cases were studied. Four analyses were combined: the search for events with 
four or more leptons plus missing energy (Section 3.1.2), the searches for tracks with large impact 
parameters (Section 3.2.2) or kinks (Section 3.3.2) and the search for tracks with an anomalous 
ionization energy loss in events with more than two tracks (Section 3.4.2). No indication of new 
physics was observed in any of the analyses. Significant overlaps exist in the signal efficiency between 
the search for four leptons plus missing energy and the large impact parameter search (up to 20 %) 
and between the large impact parameter and the kinked track search (up to 15%). The efficiencies 
in the overlap of other selections were found to be less than 0.1%. No data event was selected in 
common by two analyses. This is compatible with the shared expected background of 0.01 events 
between the search for multi-leptons plus missing energy and the large impact parameter search and 
no expected event for the other combinations. The results of the different analyses were combined 
assuming a f3/s dependence of the production cross-sections at the various centre-of-mass energies. 
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Limits on the production cross-section of neutralino pairs at ^/s = 208 GeV in the stau NLSP case 
are shown in Figure 8. In (a) the Umit vahd for ah NLSP Hfetimes is given. In most regions, cross- 
sections larger than 0.1 pb can be excluded, except for low mass differences between the neutralino 
and the stau. Figure 8 (b) shows the NLSP lifetimes at which the maximum cross-section limit is 
reached. For small mass differences between the neutralino and the stau, the limit is set by the 
short lifetime searches, whereas the searches for intermediate lifetimes set the limit in the rest of the 
parameter space. Figure 8 (c) shows the excluded cross-section for a short-lived NLSP and cross- 
sections larger than 0.1 pb can be excluded. Figure 8 (d) shows the exclusions of values larger than 
0.05 pb for a long-lived NLSP. 

In the slepton co-NLSP scenario, slightly better results are obtained because of higher efficiencies 
of the short lifetime analyses for muons and electrons compared to tau searches. As shown in 
Figure 9 (a), cross-sections larger than 0.1 pb can be excluded for all NLSP lifetimes, even for small 
mass differences between the neutralino and the slepton. Figure 9 (b) illustrates that the limits are 
set by the medium lifetime searches in most cases, while for small mass differences they are set by 
the short lifetime searches. Figures 9 (c) and (d) give the excluded cross-sections for a short-lived 
NLSP and a long-lived NLSP, with typical excluded cross-sections larger than 0.1 pb and 0.05 pb, 
respectively. 

6.1.3 Cross-section limits for smuon and selectron pair-production (stau NLSP) 

To search for pair-produced selectrons and smuons which decay via a (virtual) neutralino to the 
stau, the NLSP, four searches were combined. Topologies with a promptly decaying NLSP were 
searched for using the analysis sensitive to four or more leptons plus missing energy (Section 3.1.2). 
Candidate events for NLSPs with a medium lifetime were identified with the large impact parameter 
(Section 3.2.3) or kink search (Section 3.3.3), while possible stable NLSPs were searched for by 
selecting tracks with an anomalous ionization energy loss in events with more than two tracks 
(Section 3.4.2). For all four analyses, the number of selected data events is compatible with the 
number expected from Standard Model sources. 

Significant overlaps for the signal efficiency exist between the search for multi-leptons plus missing 
energy and the large impact parameter search (up to 26%), between the multi-lepton search and 
the kink search (up to 10 %) and between the large impact parameter and the kink search (up to 
7%). For the expected background, an overlap of 0.01 events was found between the search for 
multi-leptons plus missing energy and the large impact parameter search. In the overlap regions of 
other analyses no event is expected to be selected in common. 

The results for the different analyses were combined assuming a f5^/s dependence of the pro- 
duction cross-sections at the various centre-of-mass energies. Figure 10 (a) shows the upper limit 
at 95% C.L. on the production cross-section of smuon pairs as a function of the NLSP mass at 
y/s = 208 GeV for any lifetime of the NLSP. Cross-sections larger than 0.4 pb can be excluded, 
independent of the neutralino and NLSP mass, except for low mass differences between the smuon 
and the stau. The limit is set at short NLSP lifetimes in most cases, while for smuon masses 
below 65GeV/c^ the maximum excluded cross section is given by searches for a stable NLSP (Fig- 
ure 10(b)). Figure 10(c) shows that for a short-lived NLSP cross-sections larger than 0.4 pb can 
be excluded. For a long-lived NLSP limits on the production cross-section of smuon pairs less than 
0.1 pb are achieved (Figure 10(d)). 

Similar values of the excluded production cross-section are achieved for selectrons as shown in 
Figure 11. 
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6.1.4 Cross-section limits for chargino pair-production 

In the scenarios with a stau or a slepton being the NLSP, the charginos decay to the NLSP plus 
a neutrino. Thus the signature visible in the detector is similar to the one expected from direct 
slepton NLSP pair-production except for an additional significant amount of missing energy which 
is taken away by the neutrino. If the NLSP decays promptly or inside the detector volume the 
same searches as for direct NLSP production can be applied. These are the searches for acoplanar 
leptons (Section 3.1.1), for tracks with large impact parameters (Section 3.2.1) and for kinked tracks 
(Section 3.3.1). In case of a long-lived NLSP the analysis described in Section 3.4.2 has to be applied 
as the two NLSP tracks are no longer back-to-back. No sign of new physics was observed by any of 
the four searches nor their overlap analyses. 

For the exclusive and overlap analyses, the efficiency functions for NLSP pair-production with a 
stau or slepton NLSP were used, taking into account that the NLSP is a secondary particle. Thus, 
the efficiency function uses the average /?7 of the NLSP which is calculated from the kinematics of 
a two-body chargino decay, as described in Section 5.1.1. It was checked using a full simulation of 
chargino events that the applied interpolating functions describe the efficiencies well. 

Cross-section limits were calculated both for the case of equal branching ratios of the chargino 
to all NLSP flavours in the slepton co-NLSP scenario and 100 % branching ratio to the stau in the 
stau NLSP scenario. The data recorded at various centre-of-mass energies were combined assuming 
a /3/s dependence of the cross-sections. 

Figure 12 (a) shows the upper limit at 95 % C.L. on the production cross-section of chargino pairs 
at ^/s = 208 GeV for any lifetime of the stau NLSP. Cross-sections larger than 0.2 pb can be excluded 
for all chargino masses, independent of the NLSP mass. The limit is set at short NLSP lifetimes, 
except for chargino masses below 65GeV/c^ or close to the kinematic limit where the maximum 
excluded cross section is given by searches for a stable NLSP (Figure 12(b)). Figure 12(c) shows 
that for a short-lived NLSP cross-sections larger than 0.2 pb can be excluded. For a long-lived NLSP 
limits on the production cross-section of chargino pairs less than 0.1 pb are achieved (Figure 12 (d)). 

Similar values of the excluded production cross-sections are found for charginos in the slepton 
co-NLSP scenario as shown in Figure 13. 

6.2 The neutralino as the NLSP 

6.2.1 Cross-section limits for neutralino pair-production 

Directly produced NLSP pairs are searched for by the acoplanar photon search, which is sensitive to 
NLSP decay lengths of about 10 cm, corresponding to a lifetime of 10~^ s. If one or both neutralinos 
decay outside the detector, there is no acceptance by any analysis. 

Figure 14 gives the excluded production cross-section for neutralino pairs at a centre-of-mass 
energy of ^/s = 208 GeV. Cross-sections larger than 0.04 pb can be excluded for masses of the 
short-lived neutralino ranging from 45GeV/c^ up to the kinematic limit. 

6.2.2 Cross-section limits for slepton pair-production 

To search for pair-produced sleptons in the neutralino NLSP scenario, two analyses were applied. 
For short NLSP lifetimes with a decay length of the neutralino up to 3 m (r < 10~^ s) the search 
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for leptons plus photons and missing energy (Section 4.1.2) has sensitivity. Long NLSP hfetimes 
(r > 10^^ s) were covered by the search for acoplanar leptons (Section 3.1.1). For each flavour of 
the slepton an optimized analysis was used. The overlap of efficiencies between the analyses is of 
the order of 15 % and a Standard Model background of less than 0.01 events is expected for the 
overlap analyses. For all searches the results from the data recorded at centre-of-mass energies of 
189 GeV to 209 GeV were used. No indication for new physics in addition to the expected Standard 
Model background was observed. The datasets at different centre-of-mass energies were combined 
assuming a /3^/ s dependence of the cross-sections. 

Figures 15, 16 and 17 show the production cross-sections excluded at 95% C.L. for selectrons, 
smuons and staus, respectively. Plot (a) of each figure gives the excluded cross-section valid for any 
NLSP lifetime as a function of the slepton and the neutralino masses. Values larger than 0.2 pb can 
be excluded for all slepton flavours, independent of the neutralino mass. In plot (b) of each figure 
it is shown that the limit is set at NLSP lifetimes where the cross-over of the efficiencies of the two 
analyses takes place, shifting toward higher lifetimes with higher NLSP masses. Plots (c) and (d) of 
each figure give the excluded production cross-sections assuming a short-lived or long-lived NLSP. 
For both lifetime cases and for all fiavours cross-sections higher than 0.1 pb can be excluded. 

6.2.3 Cross-section limits for chargino pair-production 

To be sensitive to all possible signatures for pair-produced charginos in the neutralino NLSP scenario, 
a large set of analyses had to be combined. For short NLSP lifetimes, the search for photons and 
missing energy plus leptons or jets is applied (Section 4.1.2), which is split into several selections, 
optimized for low and high track multiplicities in the event. The high multiplicity selection itself 
is divided into three parts, optimized for different mass differences between the chargino and the 
NLSP. Long lifetime neutralinos in chargino candidate events are selected by the search for acoplanar 
leptons in case of leptonically decaying W bosons. For semi-leptonic and hadronic W decays the 
chargino/neutralino searches with and without an identified lepton were applied. 

For the acoplanar lepton search and for the search for photons and missing energy plus leptons 
or jets data at centre-of-mass energies ^/s = 189-209 GeV were analyzed. For the chargino searches 
with a semi-leptonic and hadronic W decay and a long-lived neutralino, data at y/s = 192-209 GeV 
were used. No indication for new physics was found for any of the analyses. Limits on the produc- 
tion cross-section for charginos were computed, with the results at different centre-of-mass energies 
combined assuming a /3/s dependence of the cross-sections and a pure decay of the chargino via a 
W boson. 

Figure 18 gives the excluded production cross-section at 95 % C.L. and for -y/s = 208 GeV for 
chargino pairs. The limits in (a) are valid for all NLSP lifetimes; cross-sections larger than 0.3 pb 
can be excluded for most chargino and NLSP masses. As can be seen in (b), this limit is set at 
NLSP lifetimes where the cross-over of the sensitivity of the lifetime analyses takes place, shifting 
toward higher lifetimes with higher NLSP mass (lower Pj of the NLSP). Figures (c) and (d) give 
the excluded production cross-section assuming a short-lived or long-lived NLSP. For short-lived 
NLSPs, cross-sections higher than 0.2 pb, for long-lived NLSPs, cross-sections higher than 0.3 pb 
can be excluded. 
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7 Interpretations in the framework of the GMSB model 



7.1 The GMSB scan database 

Interpreting the experimental results in terms of the GMSB model requires a comparison with the 
theoretical expectations within the framework of the model. Chosen here is the minimal version 
of the GMSB model with five parameters and a sign in addition to the SM parameters. The new 
parameters are the SUSY breaking scale, ^/F, the messenger scale, M, the messenger index, A^, 
the ratio of the vacuum expectation values of the two Higgs doublets, tan /3, the sign of the Higgs 
sector mixing parameter, sign(//), and the mass scale A, which determines the SUSY particle masses 
at the messenger scale. A certain point in the parameter space of the model is excluded if the 
experimental upper limit on the cross-section at ^/s = 208 GeV discussed in Section 6 is less than 
the expected cross-section a at ^/s = 208 GeV, taking into account the branching ratio BR (cjxBR^ 
for pair-produced particles decaying to the same final state). 

The parameter of the SUSY breaking scale ^/F is eliminated as experimental upper limits were 
calculated on the production cross-section of SUSY particles independent of the NLSP lifetime 
(Equation 1). In the phase space of the remaining parameters a scan was performed to calculate the 
complete mass spectrum, production cross-sections and branching ratios for different SUSY particles 
at each point considered. For this scan the framework and formulae of [4] were used and generalized 
to include a full mass treatment for all three generations. The calculations are embedded in the 
SUSYGEN generator. 



Parameter 


Scan points 


Step size 


A 


5 - 150TeV/c2 


lTeV/c2 


tan P 


1 - 50 


0.2 


M 


1.01 • A, 250TeV/c2, 10^ TeV/c^ 




N 


1, 2, 3, 4, 5 




sign(^) 


+1, -1 




Mtop 


175 GeV/c2 





Table 4: Scanned points in the GMSB parameter space. The parameter A sets the overall mass 
scale of the SUSY particles, tan/3 is the ratio of the vacuum expectation values of the two Higgs 
doublets, M is the messenger scale, the messenger index, and sign(/x) is the sign of the Higgs 
sector mixing parameter. 

The model parameters as well as the range and step size considered for them are summarized 
in Table 4. The messenger scale M is arbitrary in the minimal model, but, as the mass rrib of the 
messenger bosons is given by = M-\/l±A/M, the relation M > A has to be fulfilled in order to 
obtain a positive messenger boson mass squared. Both models with M ~ A and M ^ A are viable; 
therefore, three scenarios for the messenger scale were studied: M very close to A (M = 1.01 • A), 
M = 250TeV/c^ and M very large (lO^TeV/c^). For the integer parameter N, values up to 
five were considered. This is adequate, since perturbativity of the gauge interactions up to the 
grand unification scale, MquT; implies N ~ 150 / In ^^^7^ [^]- Thus, for a messenger mass scale 
M = lOOTeV/c^, iV < 5 is required. Both signs of the parameter n were considered. For each of the 
30 combinations of A^, M and sign(/i), a scan was performed in A and tan/3. The upper and lower 
limits for this scan were chosen according to the following considerations. For A > 150TeV/c^ the 
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Supersymmetric particles are already very heavy and cannot be produced at LEP energies. Both 
the regions with tan/3 > 50.0 and A < 5TeV/c^ are theoretically forbidden. The exact shape of the 
theoretically inaccessible region in the A — tan/3 plane depends on the other parameters. 

The scan information was saved in a database for comparison with the experimental results. In 
total, 270 such data sets were produced: 30 (with different N, M and sign(/_f)) for each centre-of- 
mass energy and for nine different energy points with y/s = 182.7, 188.7, 191.6, 195.5, 199.5, 201.6, 
205.1, 206.7 and 208.1 GeV. 

A gravitino mass of 2eV/c^ was chosen, corresponding to a SUSY breaking scale of \/F ~ 
lOOTeV/c^. This is motivated by the requirement that the branching ratio of the next-to-NLSP to 
the gravitino is small and only the NLSP decays to the gravitino. As long as this is fulfilled, the 
cross-sections and branching ratios do not depend on the gravitino mass. This makes it possible to 
decouple the issue of NLSP lifetime, which depends on VF, from the scan. Note that ^/F does not 
have a large effect on the other sparticle masses. 

7.2 Direct constraints on the NLSP mass (stau or slepton) 

Constraints on the NLSP masses were determined from the cross-section limits obtained for NLSP 
pair-production in the slepton co-NLSP scenario and stau pair-production in the stau NLSP scenario 
(Section 6.1.1). This was done by comparing the excluded cross-section at = 208 GeV with the 
production cross-section predicted by the theory. For pair-produced particles the expected cross- 
section times the branching ratio squared, a • BR^, within the model has to be known, but naturally 
this varies strongly within the model, depending on the parameter set chosen. To obtain values 
which are valid for each of the parameter sets considered, the following, conservative, minimization 
procedure was applied. 

For each of the 30 parameter combinations of M, N and sign(/i), a scan over the parameters 
A and tan/3 was performed. In the regions where the sleptons or stau are the NLSP the minimum 
(7-BR^ for each NLSP mass is determined within the parameter set considered. For selectrons strong 
variations were found because of additional positive and negative interfering t-channel production 
which contributes differently for the various parameter sets. This is in constrast to selectron searches 
in gravity mediated SUSY models in which the selectron cross-section is always enhanced by an 
interfering t-channel contribution. Finally for each NLSP mass the minimum of a ■ BR^ within these 
30 scenarios was calculated and the resulting minimal expected cross-sections were compared to the 
experimentally achieved limit. 

Figure 19 shows the mass limits for staus in the stau NLSP scenario (a), and smuons and 
selectrons in the slepton co-NLSP scenario (b,c) as a function of their lifetime. For each slepton 
flavour, the lowest mass constraints are found for very short slepton lifetimes, with the exception of 
the selectrons, where the region around 60 GeV/c^ cannot be excluded at long lifetimes. This is due 
to the effect of similar values of the ionization energy loss of heavy and light charged particles for 
a region of momenta around 65 GeV/c, thus a decreased efficiency of the searches using the dE/dx 
measurements. In the stau NLSP scenario, staus with masses below 87.4GeV/c^ are excluded at 
95% C.L. for an expected limit from the background-only hypothesis of ruf^ > 87.6GeV/c^. In the 
slepton co-NLSP scenario the directly excluded selectron, smuon and stau masses (and expected 
limits), valid for any slepton lifetime, are rngj^ > 60.1GeV/c^ (60.0 GeV/c^), m^j^ > 93.7GeV/c^ 
(93.6GeV/c2) and ruf^ > 87AGeV/c^ (88.2 GeV/c^). In the slepton co-NLSP scenario the sleptons 
are mass-degenerate; their mass difference is - by definition - less than the r, /i or e mass, so 
that decays i' — > i'i are forbidden. Therefore the highest mass limit, the limit on the smuon mass 
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of m^jj > 93.7GeV/c^ can be used to deduce an indirect common mass limit on all sleptons of 
= rrifi^ — Mr > 91.9GeV/c^ within the slepton co-NLSP scenario. 

No direct mass limits were calculated for SUSY particles heavier than the NLSP. The expected 
production cross-section for neutralinos, charginos, heavy selectrons or smuons depends strongly on 
the model parameters and is suppressed in some regions of the parameter space. Thus the mini- 
mization method gives no sensible general expectation which can be compared to the experimental 
results. 

7.3 Direct constraints on the NLSP mass (short-lived neutrahno) 

For the neutralino NLSP, the same minimization method described in the previous section (7.2) 
was applied to obtain an expected production cross-section valid at each scan point. In Figure 14 
this prediction is compared to the experimentally excluded production cross-section for neutralino 
pairs at a centre-of-mass energy of = 208 GeV. Neutralino masses less than 96.8GeV/c^ can be 
excluded at 95% C.L. for an expected limit from the background-only hypothesis of 96.3GeV/c^. 
The limit is valid for short-lived neutralinos with lifetimes up to 10^^ s only. 

In the neutralino NLSP scenario no direct mass limits were computed for charginos and sleptons. 
The production cross-sections of these particles depend strongly on the model parameters and reach 
zero in some parts of the GMSB parameter space. Thus no general minimal expectation can be 
computed to compare to the experimental results. 

7.4 Exclusions within the GMSB parameter space 

For each of the 30 sets of the GMSB parameters A^, M and sign(/i) considered in the scan, the 
exclusion in the A-tan /? plane was studied. At each point in the plane, the cross-sections for the 
various channels as well as their branching ratios are known. A point in the parameter space is 
excluded if it is kinematically accessible and the expected cross-section in at least one channel is 
higher than the experimentally derived 95% C.L. cross-section limit in this channel, taking into 
account the branching ratio(s). 

As examples, in Figure 20 the excluded regions in the A-tan j3 plane, valid for any NLSP lifetime, 
are shown for = 1 and 3, M = 1.01 • A and 250TeV/c^ and sign(/x)> 0. The shaded or hatched 
regions correspond to different search channels, indicating the relevance of the various analyses in 
the parameter space. 

In general, the stau NLSP and slepton co-NLSP scenarios become more important as number 
of messenger sets N increases. Searches for direct NLSP pair-production exclude almost completely 
the accessible region of the stau NLSP and slepton co-NLSP scenarios. Only a few additional points 
are excluded by the neutralino searches, for example regions with tan /? w 12 and A 55 TeV/c^ for 
N = 1, M = 1.01 • A in Figure 20. 

The neutralino NLSP scenario plays an important role only for N = 1 and in regions with high M 
for N = 2. In this scenario, a large fraction of the parameter space can be excluded by searches for 
slepton pair-production. Searches for chargino pair-production exclude additional points in regions 
with low tan/? and high M. In the chargino channel a complication arises due to the fact that the 
chargino has two decay modes which can lead to the same final state: 

(1) e+e- -^xtXi^ X? W+* x? W"* ^ xll^ux\^~y, 
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(2) e+e- ^ xtxl ^ i+'^ti^ ^ i+j^xl^i^xl- 



Since in the scan for the branching ratios of the different final states no information about the decay 
chain is available, it is impossible to decide if the branching ratio to the i'^v final state is due to 
leptonic W decays (1), or due to the direct decay mode (2). For the interpretation, only results from 
the decay channel (1) were used. The analyses sensitive to the leptonic final state are the search 
for acoplanar leptons in case of long neutralino lifetimes and the search for photons, missing energy 
plus leptons for short lifetimes. It was checked that these analyses have a similar efficiency for both 
chargino decay modes that lead to the leptonic final state. The branching ratio to iu was treated 
as if it was purely due to channel (1). Since the efficiency for (2) is the same, even if this branching 
ratio was totally due to decay mode (2), the cross-section limits would not be affected. 

In general, the production cross-section for charginos is larger than that for sleptons in the 
regions excluded by this channel. Nevertheless, the slepton searches exclude larger regions than 
chargino searches because in GMSB models charginos are mostly heavier than the (right-handed) 
sleptons. Thus, the charginos are kinematically limited and higher values of A can be reached by 
the slepton searches. 

7.5 Constraints on the SUSY particle mass scale A 

From the exclusions in the A-tan /? plane for each of the 30 sets of the GMSB parameters A^, M 
and sign(/_f) considered in the scan, lower limits were inferred on the SUSY particle mass scale A, 
independent of tan /?. For fixed A^, the parameter A determines the GMSB particle spectrum at the 
messenger scale, since the gaugino masses are given by 



where Aj are the gaugino fields of gauge group i and the are the GUT-scale normalized coupling 
constants of these gauge groups. Also the scalar masses at the messenger scale are determined by A: 



where the sum is over the gauge groups SU{3)c, SU{2)i and C/(l)y, and ki are constants of 0(1). 

The results, valid for all NLSP lifetimes and thus independent of the SUSY breaking scale \fF, 
are summarized in Table 5. It can be seen that the constraints on A decrease with larger A^ and 
lower limits of typically 40TeV/c^ are found for N = 1 and 15TeV/c^ for N = 5. The constraints 
depend on M, but are almost independent of sign(^). The limits are set at values of tan/3 between 
5 and 25; for higher and lower tan /? values they get larger. In conclusion, constraints on A of 
A > 40, 27, 21, 17, 15TeV/c2 were derived for A^ = 1, 2, 3, 4, 5, respectively, for all M, tan/3, 
sign(^) and all NLSP lifetimes (all values of ^/F). 

The constraints on A imply a lower limit on the neutralino mass in the neutralino NLSP scenario, 
independent of the neutralino lifetime. Masses below 53.5GeV/c^ for N=l up to 94.0GeV/c^ for 
N=5 can be excluded for all neutralino lifetimes. 



a^{M) 




39 



Lower limit at 95% C.L. on A [TeV/c^] for all NLSP lifetimes 




= 1 


N = 2 


N = 3 


N = 4 


A^ = 5 


High M {M = 10^ TeV/c2) 


/X > 


40 

43 


27 
27 


21 
21 


17 
17 


15 
15 


Medium M (M = 250TeV/c2) 


/X < 
/i > 


49 
49 


31 
34 


26 
26 


22 
22 


20 
20 


Low M (M = 1.01 • A) 


/i < 
/X > 


54 
52 


37 
37 


30 
30 


25 
25 


22 
22 



Table 5: Lower limits at 95 % C.L. on the SUSY particle mass scale A, for various sets of the GMSB 
parameters M, A^ and sign(/x), valid for all NLSP lifetimes. These constraints use only the direct 
searches for SUSY particle production. The minimum value of A for each A^ is highlighted. 

7.6 Impact of searches for the neutral Higgs boson 

The limit on the Standard Model Higgs boson reached by the four LEP experiments [33] can be 
used to place additional constraints on the GMSB parameter space. This is justified as in GMSB 
there is almost no parameter space in which the Higgs-Z-Z coupling is suppressed compared to its 
Standard Model value. Also, in GMSB the branching ratio oi H ^ bb, the most important channel 
in the SM Higgs search, shows no strong dependence on tan/3 and is suppressed by at most 10% in 
a few regions of the phase space [34]. 

In Figure 20, in addition to the regions excluded by direct sparticle searches, constraints from the 
LEP combined Higgs limit of 114.4 GeV/c^ are also shown. The Higgs constraints initially appear 
to be much stronger than those from the direct SUSY particle searches; however, for a given set of 
GMSB parameters, the theoretical uncertainty on the inferred Higgs mass is about 3 GeV/c^, and 
the Higgs mass uncertainty due to the uncertainty on Mtop is about 5 GeV/c^ [34]. The effect of 
these uncertainties is also illustrated in Figure 20, where it can be seen that the Higgs constraints can 
be rather weak when the full uncertainty is taken into account. Because of these large uncertainties, 
constraints on the GMSB parameter space from the Higgs search are not included in our quoted 
limits. 

8 Conclusions 

Searches have been performed for topologies predicted by GMSB models. All possible lifetimes of 
the NLSP which is either the lightest neutralino or a slepton, have been considered. 

No evidence for new physics has been found in the OPAL data sample collected at centre- 
of-mass energies y/s = 189-209 GeV. For the first time limits are presented on the production 
cross-section for all search topologies. The impact of the searches on the minimal GMSB model 
has been tested by performing a scan over its parameters. NLSP masses below 53.5 GeV/c^ in the 
neutralino NLSP scenario, below 87.4GeV/c^ in the stau NLSP scenario and below 91.9GeV/c^ in 
the slepton co-NLSP scenario can be excluded at 95% C.L. for all lifetimes of the NLSP. The scan 
gives constraints on the universal SUSY mass scale A from the direct searches for SUSY particle 
production of A > 40, 27, 21, 17, 15TeV/c2 for A^ = 1, 2, 3, 4, 5 for ah NLSP lifetimes. 
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Figure 1: Distributions of the event visible energy (a, b) and the number of low multiplicity jets 
consistent with having originated from leptons (c, d) in the search for multi-leptons with missing 
energy described in Section 3.1.2. The distributions are plotted with all other cuts applied. In (a) 
and (c) the points represent the selected data events from the entire 189-209 GeV data sets, while 
the shaded regions show the expectation from Standard Model backgrounds. Plots (b) and (d) show 
examples of distributions with arbitrary scale for signal events for e"'"e~ — > XiXi ~^ tttt with the 
prompt decay f — > rG at ^/s = 206.0 GeV . The solid, dashed and dotted lines correspond to (xi^t^) 
mass combinations of (50,45), (85,65) and (102,45) GeV/c^, respectively. The arrows indicate the 
cut positions. 
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OPAL, Vs = 189 - 209 GeV 
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Figure 2: Distributions of variables in the search for charged NLSPs with medium hfetime (tracks 
with large impact parameters) before a specific selection cut. The particular examples were taken 
from the search for neutrino pair-production, Section 3.2.2. The data, taken at centre-of-mass 
energies of 189-209 GeV are represented by dots, the simulated background from Standard Model 
sources by a shaded histogram. Plot (a) shows the distribution of the transverse momentum 
of the primary tracks before cut (3). Events with at least two primary tracks with a momentum 
Pt greater than l.OGeV/c are selected. The value is indicated by a hatched line and arrow. Plot 
(c) gives the invariant mass W-j^-^ of the secondary tracks assuming a decay into two charged pions 
before cut (5). The arrow gives the region which is selected by the analysis. In (b) and (d) the 
solid and dashed lines indicate the distributions with arbitrary normalization for an expected signal 
which corresponds to (x?, t^) mass combinations of (102, 100) and (50, 45) GeV/c^. 
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OPAL, Vs = 189 - 209 GeV 
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Figure 3: Distributions of variables in the search for charged NLSPs with medium hfetime (kinked 
tracks) before a specific selection cut. The particular examples were taken from the search for direct 
NLSP pair-production, Section 3.3.1. The data, taken at centre-of-mass energies of 189-209 GeV are 
represented by dots, the simulated background from Standard Model sources by a shaded histogram. 
In (b, d) the solid, dashed and dotted lines indicate the distributions expected from a pair-produced 
stau with masses of 65, 80 and 95GeV/c^. Plots (a, b) show the distribution of the transverse 
momentum of the primary tracks before cut (6). Plots (c, d) give the invariant mass Wqo of the 
primary track assuming a decay into two massless particles before cut (7). The arrows indicate the 
regions which are selected by the analysis. 
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OPAL, Vs = 130 - 208 GeV 
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Figure 4: Distributions of variables in the search for long-hved charged NLSPs before selection cuts. 
Plot (a) shows the scaled component Pz of the event momentum P = |P| = | X^Ptrackl along the beam 
axis before cut (5), Section 3.4.2. The data, taken at centre-of-mass energies of 189-209 GeV are 
represented by dots, the simulated background from Standard Model sources by the solid line. The 
dashed line indicates the expected distribution for staus with a mass of 90 GeV/c^, pair-produced at 
a centre-of-mass energy of 208 GeV. In Plot (b) the distribution of candidate events before cut (9) 
as a function of the track momentum P and the specific ionization energy loss dE/dx is given. In 
addition the expected energy loss for electrons (e), muons (/i), pions (tt), kaons {K)^ protons (p) 
and long-lived staus (fi with mf^ = 45 and 90GeV/c^) is indicated. In both plots the arrows give 
the regions which are selected by the analysis cuts. 
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OPAL, Vs = 206 GeV 
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Figure 5: Distributions of variables in the search for slepton and charginos with very short to medium 
lifetime (Section 4.1.2), selection (A). Plots (a, b) show the distribution of the visible energy in the 
event as a fraction of the centre-of-mass energy, plots (c, d) the energies of the most energetic 
photon in the event. In (a) and (c) the data recorded at ^/s = 206 GeV are represented by dots, the 
expected Standard Model background by the shaded histogram. In (b) and (d) the distributions 
of a possible signal of stau pair-production in the neutralino NLSP scenarios are given. The solid, 
dashed and dotted lines correspond to (fj^, Xi) mass combinations of (102, 51), (80, 20) and (50, 
45) GeV/c^. The arrows indicate the cut positions. 
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Figure 6: Efficiencies for stau pair-production at ^/s = 208 GeV, as a function of the lifetime 
for fixed stau masses. The symbols represent the efficiencies for ten simulated lifetimes and the 
curves show the interpolating efficiency functions for the different searches: the exclusive search for 
promptly decaying staus (dashed), the exclusive search for large impact parameters (long dashed- 
dotted), the exclusive search for kinks (dotted) and the exclusive search for stable staus (dashed- 
dotted). The overlap efficiencies between these searches are shown as filled histograms. The sum of 
all efficiencies, i.e. the four exclusive selection efficiencies and the four overlap efficiencies, is shown 
by the squares and the corresponding solid line. 
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Figure 7: Contours of the 95% C.L. upper limits on the production cross-sections (T95 for slepton 
pair-production at ^/s = 208 GeV as a function of the slepton mass and lifetime riifc. Shown are 
cross-section limits for (a) staus in the stau NLSP scenario, (b) staus, (c) smuons and (d) selectrons 
in the slepton co-NLSP scenario. The shadings correspond to different ranges of the upper limit on 
the cross-section, as indicated by the scale on the right side . 
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Figure 8: Upper limits at 95 % C.L. on the production cross-section for neutralino pairs in the stau 
NLSP scenario at ^/s = 208 GeV. The data analyzed, taken at different centre-of-mass energies, 
were combined assuming a, (3/ s dependence of the cross-sections. In (a) the lifetime independent 
exclusions are shown as a function of the neutralino and NLSP mass. The lifetime at which the limit 
is set is given in (b). Figures (c) and (d) give the excluded production cross-section for neutralinos 
assuming a very short-lived or a stable NLSP. 
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Figure 9: Upper limits at 95 % C.L. on the production cross-section for neutralino pairs in the slepton 
co-NLSP scenario at ^/s = 208 GeV. The data analyzed, taken at different centre-of-mass energies, 
were combined assuming a, (3/ s dependence of the cross-sections. In (a) the hfetime independent 
exclusions are shown as a function of the neutrahno and NLSP mass. The hfetime at which the hmit 
is set is given in (b). Figures (c) and (d) give the excluded production cross-section for neutralinos 
assuming a very short-lived or a stable NLSP. 



51 



OPAL, eV |Xr|Xr (Tj NLSP) 



a) all NLSP lifetimes 




100 

m{\i^) (GeV/c) 



b) NLSP lifetime at max. a, 



log(x,=f,/s) 



4i 100 

> 

O 90 

'S" 80 
B 

70 
60 
50 



95 






80 



100 

m(|Ig) (GeY/c) 



I 



life' 

-6 

-7 

-8 

-9 

-10 

-11 

-12 




d) NLSP lifetime > 10 



100 

mdlg) (GeV/c^) 




100 

mdl^) (GeV/c^) 



Figure 10: Upper limits at 95% C.L. on the production cross-section for smuon pairs in the stau 
NLSP scenario at ^/i = 208 GeV. The data analyzed, taken at centre-of-mass energies of y/s = 189- 
209 GeV, were combined assuming a f3^/s dependence of the cross-section. In (a) the lifetime 
independent exclusions are shown as a function of the smuon and NLSP mass. The lifetime at 
which the limit is set is given in (b). Figures (c) and (d) give the excluded production cross-section 
for smuons assuming a very short-lived or a stable NLSP. 
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Figure 11: Upper limits at 95% C.L. on the production cross-section for selectron pairs in the stau 
NLSP scenario at ^/s = 208 GeV. The data analyzed, taken at centre-of-mass energies of ^/s = 189- 
209 GeV, were combined assuming a dependence of the cross-section. In (a) the lifetime 

independent exclusions are shown as a function of the selectron and NLSP mass. The lifetime at 
which the limit is set is given in (b). Figures (c) and (d) give the excluded production cross-section 
for selectrons assuming a very short-lived or a stable NLSP. 
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Figure 12: Upper limits at 95% C.L. on the production cross-section for chargino pairs in the stau 
NLSP scenario at y/s = 208 GeV. The data analyzed, taken at different centre-of-mass energies, 
were combined assuming a (3/s dependence of the cross-sections. In (a) the lifetime independent 
exclusions are shown as a function of the chargino and the NLSP mass. The lifetime at which 
the limit is set is given in (b). Figures (c) and (d) give the excluded production cross-section for 
charginos assuming a very short-lived or a stable stau. 
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Figure 13: Upper limits at 95 % C.L. on the production cross-section for chargino pairs in the slepton 
co-NLSP scenario at ^/s = 208 GeV. The data analyzed, taken at different centre-of-mass energies, 
were combined assuming a, (3/ s dependence of the cross-sections. In (a) the lifetime independent 
exclusions are shown as a function of the chargino and the NLSP mass. The lifetime at which 
the limit is set is given in (b). Figures (c) and (d) give the excluded production cross-section for 
charginos assuming a very short-lived or a stable NLSP. 
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Figure 14: The excluded production cross-section at ^/s = 208 GeV for neutralino pairs in the 
neutrahno NLSP case is given by the grey region. The dotted hne represents the expected hmit. 
Both hmits are given at 95 % C.L. The sohd hne shows the minimum cross-section for neutrahno pairs 
predicted by the theory. Within the model, neutralino masses between 45GeV/c^ and 96.8GeV/c^ 
(96.3GeV/c^ expected) can be excluded. The result is valid for promptly decaying NLSPs only. 
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Figure 15: Contours of the cross-section limits at 95% C.L. and y/s = 208 GeV for selectron pair- 
production in the neutrahno NLSP scenario. The cross-section Umits in (a) are vahd for any neu- 
trahno hfetime. Plot (b) shows the lifetimes at which the limit is set. The 95% C.L. cross-section 
limits for prompt neutralino decays, where only the search for photons, missing energy plus lep- 
tons/jets contributes, are plotted in (c). In (d) the 95% C.L. cross-section limits are shown for very 
long-lived neutralinos, where only the search for acoplanar leptons contributes. 
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Figure 16: Contours of the cross-section limits at 95% C.L. and y/s = 208 GeV for smuon pair- 
production in the neutrahno NLSP scenario. The cross-section Umits in (a) are vahd for any neu- 
trahno hfetime. Plot (b) shows the lifetimes at which the limit is set. The 95% C.L. cross-section 
limits for prompt neutralino decays, where only the search for photons, missing energy plus lep- 
tons/jets contributes, are plotted in (c). In (d) the 95% C.L. cross-section limits are shown for very 
long-lived neutralinos, where only the search for acoplanar leptons contributes. 
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Figure 17: Contours of the cross-section limits at 95% C.L. and y/s = 208 GeV for stau pair- 
production in the neutrahno NLSP scenario. The cross-section Umits in (a) are vahd for any neu- 
traUno hfetime. Plot (b) shows the lifetimes at which the limit is set. The 95% C.L. cross-section 
limits for prompt neutralino decays, where only the search for photons, missing energy plus lep- 
tons/jets contributes, are plotted in (c). In (d) the 95% C.L. cross-section limits are shown for very 
long-lived neutralinos, where only the search for acoplanar leptons contributes. 
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Figure 18: Contours of the cross-section limits at 95% C.L. and y/s = 208 GeV for chargino pair- 
production in the neutralino NLSP scenario assuming a 100% branching faction for — > W^*x''- 
The cross-section hmits in (a) are vahd for any neutrahno hfetime. Plot (b) shows the NLSP lifetimes 
at which the maximum excluded cross-section is found. The results for prompt neutrahno decays, 
where mainly the searches for photons and missing energy plus leptons/jets contribute, are plotted 
in (c). In (d) the limits are shown for very long-lived neutralinos. 
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Figure 19: The observed lower mass limits for pair-produced staus in the stau NLSP (a) and 
smuons (b), selectrons (c) in the slepton co-NLSP scenario as a function of the particle lifetime 
using the direct search. For staus in the slepton co-NLSP scenario the observed and expected 
lower limit are identical to the limits of the stau in the stau NLSP scenario. The mass limits are 
valid for a messenger index N< 5. For the stau NLSP and slepton co-NLSP scenarios, the NLSP 
mass limits are set by the stau mass limit (m-NLSP > 87.4GeV/c^ (a)) and by the smuon mass limit 
(jTiNLSP > 93.7GeV/c^ (b)), respectively. 
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Figure 20: Examples of excluded regions in the A — tan /? plane excluded by pair-production searches 
for different particles, with sign(/i) > and valid for any NLSP lifetime for four different sets of 
parameters, = 1 or 3 and M = 1.01 • A or 250TeV/c2. 
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